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1.  Analysis  of  potential  schemes  of  nonequilibrium  MHD  generator  in  scramjet 
inlet  and  consideration  of  theirs  realizability  will  be  done.  Possible  configurations 
of  magnetic  field,  location  and  configuration  of  electrodes  will  be  discussed.  More 
preferable  configuration  will  be  chosen 

The  first  stage  of  the  project  is  devoted  to  analysis  of  potential  schemes  of  the 
nonequilibrium  MHD  generator  in  the  scramjet  inlet,  possible  configurations  of  a  magnetic  field, 
location  and  configuration  of  electrodes.  The  main  objective  for  the  nonequilibrium  MHD 
generator  in  the  scramjet  inlet  is  improvement  of  the  inlet  characteristics  at  off-design 
conditions.  It  is  well  known,  that  at  flight  Mach  numbers  Moo  exceeding  design  Mach  number 
Md,  the  MHD  generator  can  be  used  for  modification  of  the  flow  field  in  the  inlet  and  making  it 
quite  similar  to  the  one  at  design  conditions.  At  flight  Mach  numbers  below  the  design  Mach 
number  the  MHD  generator  allows  one  to  increase  the  air  capture  in  the  inlet. 

In  order  to  ensure  significant  MHD  influence  on  a  flow  in  the  inlet,  it  is  necessary  to 
create  a  magnetic  field  and  provide  high  enough  level  of  an  electric  conductivity  of  the  flow. 
The  static  temperature  of  the  flow  in  the  inlet  is  insufficient  to  provide  equilibrium  conductivity, 
thus  it  is  necessary  to  ionize  the  flow  using  methods  for  producing  nonequilibrium  conductivity. 
In  our  project  the  electron  beam  is  considered  as  ionizer.  As  a  base  for  the  inlet  geometry,  the 
two-shock  inlet  with  the  total  turning  angle  of  15°  and  the  Mach  design  number  Md=10  is 
chosen.  The  scheme  of  the  inlet  is  shown  in  Fig.l.  The  total  turning  angle  in  the  inlet  is 
0n=0i+02=1  5°,  where  0i=6.5°  and  02=8.5°,  the  height  of  the  throat  is  Fth=0.12m.  Flow  field  in 
the  MHD  controlled  inlet  will  be  calculated  in  2D  Euler  approach  for  selected  configurations  of 
magnetic  field  and  ionized  region.  Computational  procedure  is  based  on  the  explicit  high- 
resolution  shock-capturing  Godunov- type  scheme. 

A  source  of  magnetic  field  in  the  project  is  electromagnet.  The  top  values  of  a  magnetic 
flux  density  can  be  obtained  if  we  use  the  superconducting  electromagnet.  It  seems  more  realistic 
to  locate  the  electromagnet  inside  the  body  of  the  scramjet  inlet,  because  placement  of  any  part 
of  the  magnet  in  a  flow  outside  the  body  can  increase  the  drag  of  the  flight  vehicle.  Feasible 
location  points  of  electromagnets  in  the  scramjet  inlet  were  analyzed  at  the  first  stage  of  the 
project.  A  version  of  the  electromagnet  coil  location  enabling  realization  of  the  MHD  influence 
on  the  flow  in  the  inlet  is  presented  in  Fig.2.  In  this  case,  the  magnetic  field  is  produced  by  a 
single  coil  which  has  a  profile  adapting  both  to  design  of  the  flight  vehicle  and  requirements  to 
the  MHD  generator  which  controls  the  flow  in  the  inlet.  In  the  figure  the  coil  is  not  closed.  Its 
prolongation  can  engage  the  isolator  which  is  positioned  upstream  of  the  combustion  chamber, 
the  combustion  chamber  itself  and  the  part  of  the  scramjet  downstream  of  it.  The  magnetic  field 
in  the  above  subsystems  of  the  scramjet  can,  in  principle,  be  used  to  control  the  internal  flow  in 
the  scramjet  in  order  to  improve  its  performance.  The  scheme  of  the  scramjet  with  MHD  bypass 
can  be  used  for  the  purpose.  The  coil  in  the  plane  of  its  cross-section  basically  can  havejust_any 
profile.  In  this  project  we  restrict  our  analysis  to  only  two  profiles  consideration,  namely 
rectangular  and  circular  ones.  Feasible  location  points  of  electrodes  are  shown  schematically  in 
Fig.2.  The  electrodes  denoted  in  the  figure  by  number  1  and  number  2  can  be  used  in  the  MHD 
controlled  inlet  both  jointly  and  separately. 

To  calculate  the  magnetic  field  in  the  region  of  MHD  interaction  we  use  the  non- 
inductive  approach,  typical  for  MHD  generator  conditions.  This  approach  corresponds  to 
magnetic  Reynolds  number  Rem«l.  In  this  case,  a  magnetic  field  induced  by  current  in  the 
plasma  is  negligible  by  comparison  with  the  magnetic  field  created  by  the  current  in  the 
electromagnet  coil.  The  magnetic  field  distribution  is  calculated  by  means  of  Maxwell's 
equations.  The  steady  state  equations  are  considered  for  which:  rotB  =  /u\c,  divB  =  0 ,  where  B 

is  the  magnetic  induction  vector,  jc  is  the  current  density  in  the  electromagnet  coil,  p  is  the 
magnetic  permeability  of  a  medium.  In  calculations  we  use  p=po=47iTO_7  N/A2.  To  determine 
the  magnetic  field  we  introduce  the  vector  potential  A,  which  is  connected  with  the  magnetic 
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field  by  the  ratio  B  =  rot  A .  Equations  for  Cartesian  components  of  the  vector  potential  have  a 
form  of  Poisson's  equation: 

AAj=Mjci,  i  =  x,y,z  (1) 

If  we  know  the  current  density  in  the  electromagnet  coil,  for  example,  function  jcz(x,y),  then 
solution  of  equations  (1),  correspondingly  for  Az,  can  be  obtained  in  the  form  [1]: 

■4  00  00  ■* 

A(x’y)  =  y-  J  f  I  2  2  d^drj  (2) 

171  -----  v(x-£)  +{y~ri) 

Equations  for  other  components  of  the  vector  potential  have  the  same  form. 

The  magnetic  flux  density  can  be  determined  from  the  vector  potential  by  using  the  well 
known  equations: 

B  = 

dy  dz 


y  dz  dx 


(3) 


B  _  dAy  dA* 
dx  dy 

It  is  shown  in  Ref.  [2]  that  a  positive  effect  of  MHD  control  in  the  MHD  controlled  inlet 
at  flight  Mach  numbers  greater  than  the  design  Mach  number  (Moo>Md)  is  achievable  in  a  wide 
range  of  variation  of  such  characteristics  as  the  magnetic  field  orientation  and  a  configuration  of 
the  region  of  MHD  interaction.  On  the  other  hand,  in  the  case  of  Moo<Md  the  effect  of  MHD 
control  depends  significantly  on  the  magnetic  field  orientation.  That  is  why  we  try  to  satisfy  the 
requirements  for  the  MHD  control  in  the  inlet  at  off-design  conditions  with  M.x<Mt|  in  choosing 
the  magnetic  field  configuration.  The  more  significant  positive  effect  of  MHD  control  in  the  inlet 
at  Moo<Md  consists  in  increasing  the  air  capture.  In  order  to  effectively  increase  the  air  capture  in 
the  MHD  controlled  inlet,  it  is  necessary  to  organize  the  MHD  influence  on  the  flow  in  such  a 
way  that  the  Lorentz  force  projection  directed  to  the  inlet  body  receives  a  maximal  value  [3].  The 
Lorentz  force  is  defined  by  the  ratio  f  =  j  x  B  .  The  current  density  in  the  region  of  MHD 


interaction  is  determined  by  the  Generalized  Ohm's  law  j  +  /./, ( j xB)  =  cr(E  +  vxB),  where  //e 


is  the  electron  mobility,  E  is  the  electric  field  intensity. 

The  scramjet  inlet  lock-on  in  the  Cartesian  coordinates  corresponds  to  one  shown  in 
Figs.  1,2.  Two-dimensional  approach  in  which  the  flow  parameters  along  the  OZ  axis  are 
invariable  is  considered.  The  nonequilibrium  MHD  generator  configuration  must  provide  high 
value  of  the  Lorentz  force  directed  to  the  inlet  body.  So  it  makes  some  requirements  to 
configuration  of  the  magnetic  field  and  the  location  of  the  region  of  MHD  interaction.  Besides, 
the  absolute  value  of  the  Lorentz  force  projected  on  the  Z  axis  \fz\  needs  to  be  minimized.  It 

follows  from  the  analysis  made  that  the  Faraday  MHD  generator  with  sectioned  electrodes  meets 
the  requirements  substantially.  In  case  of  ideally  sectioned  Faraday  MHD  generator  fz=0.  As  for 
space  position  of  electrodes,  the  magnetic  field  nonuniformity  allows  one  to  locate  the  electrodes 
both  on  the  surface  of  the  flight  vehicle  and  in  the  plane  located  at  some  angle  to  it. 
Schematically  location  of  the  electromagnet  coil  and  MHD  electrodes  and  also  configuration  of 
the  magnetic  field  and  the  MHD  current  are  shown  in  Figs. 3,4.  Fig.3  corresponds  to  the  region 
with  electrodes  which  is  marked  by  number  2  in  Fig.2b.  Fig.4  corresponds  to  the  region  with 
electrodes  which  is  marked  by  number  1  in  Fig.2b. 

A  part  of  the  MHD  controlled  inlet  with  two  types  of  conductor  producing  a  magnetic 
field  is  shown  schematically  in  Fig.5.  Here  the  conductors  with  circular  and  rectangular  cross- 
sections  which  carry  the  current  Ic  along  the  axis  OZ  are  shown.  Parameters  which  determine  the 
value  and  the  space  distribution  of  the  magnetic  field  are  shown  in  Figs. 5a, b.  Here  we  are  listing 
the  following  parameters:  position  (xo,yo)  of  a  conductor  center,  magnitude  of  the  current  Ic, 
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width  h  and  height  d  for  the  rectangular  conductor  and  also  angle  0C  for  orientation  of  the 
conductor  plane  relatively  to  the  x  axis.  In  case  of  a  conductor  with  a  circular  cross  section  a 
magnetic  field  outside  the  conductor  depends  upon  the  distance  from  the  point  (xo,yo)  and  the 
value  of  Ic  and  doesn’t  depend  on  the  radius  of  the  conductor.  It  is  shown  in  calculations,  that  in 
case  of  a  rectangular  conductor  with  realistic  proportions  in  the  MHD  controlled  inlet  (d«h)  the 
magnetic  field  distribution  practically  doesn’t  depend  upon  the  d  value.  The  magnetic  field 
distributions  in  the  MHD  controlled  inlet  for  conductors  with  rectangular  and  circular  cross- 
sections  are  compared  in  Fig.6.  It  is  easy  to  see  that  near  the  inlet  surface  the  magnetic  fields  for 
the  cases  significantly  differ.  When  moving  away  from  the  surface,  the  difference  between  the 
magnetic  fields  is  less. 

We  consider  the  electron  beam  as  an  ionizer  which  produces  the  nonequilibrium 
conductivity  of  the  flow  in  the  scramjet  inlet.  Estimations  show  that  in  typical  conditions  of  the 

mv 

MHD  controlled  inlet  a  magnitude  of  the  Larmor  radius  rL  =  —  for  electrons  of  the  e-beam  is 

eB 

significantly  less  than  the  path  length  of  the  electrons  between  elastic  collisions,  here  m  is  the 
electron  mass,  v  is  a  velocity  of  the  electron,  e  is  the  electron  charge.  Therefore,  the  electrons  in 
the  e-beam  propagate  in  the  flow  practically  along  the  magnetic  field  line.  Depth  of  propagation 
of  electrons  in  the  flow  depends  upon  the  flow  density  and  the  initial  energy  of  the  electrons.  So, 
a  location  and  proportions  of  the  region  of  the  nonequilibrium  conductivity  depend  on  the 
magnetic  field  configuration,  the  space  location  of  the  e-beam  emitted  into  the  flow,  the  e-beam 
characteristics  and  the  flow  characteristics.  Possible  space  location  of  e-beam  in  the  MHD 
controlled  inlet  is  shown  in  Fig.7.  This  situation  corresponds  to  the  magnetic  field  produced  by 
the  conductor  with  a  circular  cross-section  centered  at  the  point  (xo,yo).  In  this  case  the  electrons 
in  the  e-beam  move  along  the  concentric  circles.  The  geometrical  locus  of  the  region  of  the 
nonequilibrium  ionization  can  be  characterized  by  the  center  position  (xo,yo),  magnitudes  of 
external  Rext  and  internal  Rint  circles  and  the  angle  0  which  determines  the  electrons  ranges  in  the 
flow.  In  order  to  determine  the  region  of  the  nonequilibrium  ionization  in  general  with  an 
arbitrary  configuration  of  the  magnetic  field,  it  is  necessary  to  solve  a  differential  equation  for 
trajectories  of  electrons  in  the  e-beam  traversing  the  inlet  surface  in  any  given  area.  We  consider 
configuration  in  which  the  magnetic  field  is  produced  by  a  current  or  a  system  of  currents 
directed  along  the  OZ  axis.  One  of  the  examples  of  such  configuration  is  shown  in  Fig.4.  As  in 
this  case  the  magnetic  field  has  only  two  components  Bx(x,y j  and  By(x,y),  the  differential 
equation  for  a  trajectory  of  electrons  from  e-beam  which  are  passing  through  a  point  (x\,y\)  has 
the  form: 

dx  Bx 

with  the  boundary  condition  y(x=x\)=y\.  The  curves  passing  through  the  points  bounding  the 
area  of  the  e-beam  emission  from  the  inlet  surface  determine  the  bounds  of  the  region  of 
nonequilibrium  ionization.  Depth  of  the  e-beam  penetration  is  considered  at  this  stage  as  a 
parameter  which  can  be  varied  in  calculations.  Power  density  q,  deposited  by  the  e-beam  in  the 
region  of  ionization  is  supposed  to  be  constant  in  the  region.  The  ionization  degree  in  the  flow 
and  its  conductivity  are  calculated  in  using  approximating  functions  from  Ref.  [4],  At  the 
following  stages  of  the  project  a  method  of  calculation  of  the  region  of  the  nonequilibrium 
ionization  will  be  developed  which  will  provide  for  more  adequate  description  of  the  real 
conditions  in  the  MHD  controlled  inlet. 

In  order  to  define  what  profile  of  a  conductor  producing  a  magnetic  field  in  the  MHD 
controlled  inlet  is  preferable,  computations  of  flow  fields  in  the  inlet  with  the  configuration 
shown  in  Fig.4  were  made.  Flow  regime  with  Moo<Md  was  considered.  As  a  more  significant 
characteristic,  the  air  mass  flow  rate  in  the  inlet  was  considered  in  analysis.  The  position  of  the 
magnetic  field  source  (xo,yo)  and  also  the  values  h  and  9C  (for  rectangular  conductor )  were  varied 
in  the  computations.  Bounds  of  the  e-beam  emission  in  the  flow  were  varied  in  the  computations 
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too.  As  a  result  of  computations  it  is_shown  that  the  MHD  control  allows  one  to  increase  the  air 
capture  in  the  inlet  in  using  both  circular  and  rectangular  profiles  for  the  conductor. 
Nevertheless,  usage  of  the  conductor  with  a  circular  cross-section  presents  a  small  advantage  in 
the  air  capture.  This  is  demonstrated  in  Fig.8,  here  cph  is  the  air  mass  flow  rate  in  the  MHD 
controlled  inlet  with  a  current  geometry  corresponding  to  Fig.5b,  (pR  is  the  air  mass  flow  rate  in 
the  MHD  controlled  inlet  with  a  current  geometry  corresponding  to  Fig.5a. 

We  have  also  investigated  characteristics  of  the  MHD  controlled  inlet  in  which  a 
magnetic  field  is  produced  by  two  antiparallel  currents  as  is  shown  in  Fig.9.  Computations  show 
that  at  such  configuration  of  a  magnetic  field  the  MHD  control  of  a  flow  in  the  MHD  controlled 
inlet  is  more  effective  than  at  configurations  of  a  magnetic  field  considered  previously.  For 
example,  in  the  configuration  shown  in  Fig.7  at  flight  Mach  number  Moo=8  and  IC=107A  a 
relative  increase  of  the  air  capture  due  to  MHD  control  doesn’t  exceed  2.5%.  At  the  same  time, 
in  the  configuration  shown  in  Fig.9  the  relative  increment  grows  up  to  10%.  Fig.  10  demonstrates 
the  flow  field  in  the  MHD  controlled  inlet  calculated  for  configuration  shown  in  Fig.9.  One  can 
see  that  a  streamline  entering  the  region  of  MHD  interaction  is  deflected  to  the  body.  This  effect 
is  more  obvious  in  the  fragment  of  the  flow  field  presented  in  Fig.  10b.  It  is  this  kind  of 
transformation  of  a  flow  that  increases  the  air  capture  in  the  MHD  controlled  inlet. 

Characteristics  of  the  MHD  controlled  inlet  computed  at  this  stage  will  be  presented  in 
the  annual  report  for  the  set  of  locations  of  the  magnetic  field  sources,  different  values  of  the 
current  Ic,  at  various  characteristics  of  MHD  generator  and  diverse  locations  of  an  area  of  the  e- 
beam  input  into  the  flow. 


a)  side  view 
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Electrodes 


b)  bottom  view 

Fig.2  Sketch  for  possible  configuration  of  an  electromagnet  coil  adapting  to  the  construction  of 
the  scramjet  inlet 
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b) 

Fig.3  Sketch  for  configuration  of  a  magnetic  field  and  MHD  currents  corresponding  to  location 
of  electrodes  marked  by  number  2  in  Fig.2b. 


Fig.4(  Sketch  for)  configuration  of  a  magnetic  field  and  MHD  currents  corresponding  to  location 
of  electrodes  marked  by  number  1  in  Fig.2b 
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a)  b) 

Fig.  5  Part  of  the  MHD  controlled  inlet  with  a  current  which  produces  a  magnetic  field 


Fig.6  Contours  of  the  magnetic  flux  density  for  two  various  currents.  Black  lines  correspond  to 
the  circular  conductor,  red  lines  correspond  to  the  rectangular  current  with  h=0.5m,  0C=8°; 
Ic=5-106A,  xo=0.4m,yo=4.1m. 
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Fig.7  Sketch  for  MHD  controlled  inlet  with  circular  magnetic  field. 


Fig.8  Relative  air  mass  flow  rate  in  the  MHD  controlled  inlet. 
Moo=6,  Md=10,  0C=9°,  x0=4.0m,  y0=0.4m,  IC=107A 
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field. 


a) 


X,  m 
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1.1 


b) 

Fig.  10  Density  contours  in  the  MHD  controlled  inlet  (a)  and  its  fragment  (b)  at  Mr/,=8,  Md=10, 
xa=4.0m,  ya=0.4m;  Xb=5.35m,  yb=0.8m,  IC=107A 


11 


2.  Kinetic  scheme  for  description  of  nonequilibrium  plasma  creation  in  channel  of 
mhd  generator  by  e-beam  ionizer  will  be  developed.  Influence  of  induced  electric 
field  in  mhd  channel  on  electron  concentration  sustained  by  e-beam  ionizer  will  be 
investigated,  electron  concentration  will  be  calculated  for  steady  and  pulse  periodic 
e-beams  both  for  one  beam  and  for  set  of  spaced  beams,  power  cost  to  creation  of 
any  given  electron  concentration  in  typical  for  mhd-controlled  inlet  conditions  will 
be  estimated 


The  second  stage  of  the  project  is  devoted  to  development  of  kinetic  scheme  for 
description  of  nonequilibrium  plasma  creation  in  channel  of  MHD  generator  by  e-beam  ionizer. 
Concentrations  of  charged  particles:  electrons,  positive  and  negative  ions  are  calculated  for 
steady  and  for  pulse  periodic  e-beams  both  for  homogeneous  ionizer  and  for  set  of  spaced 
beams.  Dependencies  of  electron  concentration  are  presented  in  function  of  power  density  spent 
on  ionization  in  typical  for  MHD-controlled  inlet  conditions.  It  allows  one  to  estimate  power 
costs  for  creation  of  any  given  electron  concentration  in  the  MHD-controlled  inlet.  Influence  of 
induced  electric  field  in  MHD  channel  on  electron  concentration  sustained  by  e-beam  ionizer  is 
investigated  too.  The  simple  approach  for  calculation  of  electron  temperature  in  conditions  of 
MHD  controlled  inlet  is  proposed. 

The  two-shock  inlet  with  the  total  turning  angle  0N=01+02=15°,  where  01=6.5°  and 
02=8.5°  was  considered  in  calculations.  In  the  MHD  controlled  inlet  there  are  three  different 
regions  of  a  flow  in  which  efficiency  of  the  e-beam  ionizer  will  be  different  too.  They  are 
denoted  in  Fig.l  by  next  symbols:  FS  is  the  free  stream,  SI  is  the  flow  after  the  first  oblique 
shock  and  S2  is  the  flow  after  the  second  oblique  shock.  Flow  parameters  in  these  regions  are 
calculated  numerically  in  using  2D  Euler  approach.  Gas  concentration  N  in  the  free  stream  flow 
(the  region  FS)  depending  on  the  flight  Mach  number  Moo  is  presented  in  Fig.2  for  two  values  of 
the  free-stream  dynamic  pressure  qco. 

In  this  stage  we  consider  approach  in  which  plasma  concentration  is  changed  only  along 
the  direction  of  the  flow  velocity  (x  in  our  designations).  Transversal  modification  of  plasma 
concentration  due  to  the  e-beam  energy  degradation  will  be  considered  in  the  following  stages. 
To  calculate  non-equilibrium  conductivity  of  a  flow  in  channel  of  MHD  generator  with  e-beam 
ionizer  we  consider  model  of  a  low  ionized  plasma  consisting  of  neutral  molecules,  electrons, 
negative  and  positive  ions  with  corresponding  number  densities:  N(X  N  ,  ne,  n+.  In  report  [1]  it 

is  shown  that  results  obtained  in  such  type  of  model  are  in  good  accordance  with  ones  obtained 
in  the  model  taking  into  account  more  than  40  plasma  components  and  more  than  230  reactions 
of  the  plasma  components.  Moreover,  according  to  Ref.  [2]  predictions  of  the  model  are  adjusted 
with  experiment.  Set  of  kinetic  equations  for  concentrations  of  the  plasma  components  is  the 
next: 


r)y\  C)Yl 

-^-  +  v-^  =  I(t,X)  +  kd-N02-n_- kaN]h ne  -  f3ein+ne 

dn,  dn^  T,  .  _ 

+  V~^T  =  I(t’  X)  “  Ptin+n-  ~  PeWe 

dn_  dn_  ,  _r2  ,  Ar  0 

“  +  =  kaN(Xne  ~  KN0 ,  •  n-  ~  Pun-n- 


dt  8x 


(1) 
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where  v  is  the  flow  velocity,  I(t,  x )  is  e-beam  induced  ionization  rate,  N0i  is  concentration  of 


oxygen  molecules  in  air,  ka  and  kd  are  the  rate  constants  of  attachment  and  detachment  of 
electrons  respectively,  (3  and  fin  are  the  rate  constants  for  electron-ion  and  ion-ion  recombination 
respectively.  These  constants  are  functions  of  gas  temperature  T  and  electron  temperature  Te.  In 
calculations  we  use  dependencies  from  papers  [3,  4],  In  kinetic  scheme  (1)  the  processes  of 
attachment  and  detachment  of  electrons  to  nitrogen  molecules  are  not  accounted  because  in 
considered  conditions  their  rate  constants  can  be  neglected  in  comparing  with  ones  for 
corresponding  processes  with  oxygen  molecules  participated  in.  Ionization  rate  /  is  determined  in 
terms  of  e-beam  characteristics  by  the  relation: 


W; 


u/w.) 

w. 


(2) 


where  jb  is  the  e-beam  current  density,  Eb  is  the  energy  of  electron,  Y  is  the  electron  stopping 
power,  p  is  the  gas  density,  Wt  is  the  ionization  cost.  Dimensionless  function  / (t,  x)  determines 
spatial  and  temporal  distribution  of  the  power  deposition  into  the  flow.  Value  qion  in  Eq.(2) 
determines  power  lost  by  electron  beam  in  the  flow  in  unit  of  time  in  unit  of  volume.  The  value 
qion  is  used  in  calculations  as  a  parameter  which  characterizes  the  ionizer.  In  the  first  phase  of 
this  stage  of  investigations  we  consider  approach  Te=T  which  is  correct  in  absence  of  magnetic 
field. 

Concentrations  of  charged  particles  were  calculated  for  different  values  of  qion,  different 
flow  parameters  and  different  functions /(t,  x) .  Firstly  the  steady  state  approach  was 

investigated  for  which  dne  +  _/dt  =  Oand  correspondingly  f(t,x)  =  f(x) .  Fig.3  shows  electron 

,  .  [0,  x  <  0 

concentration  calculated  in  considering  the  e-beam  occupying  the  half-space:  /  (x)  =  j  ^  ^  . 

One  can  see  that  in  this  case  electron  concentration  increases  while  increasing  the  x  value  and 
tends  to  a  limiting  magnitude.  These  limiting  magnitudes  of  electron  concentration  are  shown  in 
Fig.4  in  function  of  the  power  density  spent  on  the  flow  ionization  for  various  values  of  the  flight 
Mach  number.  In  order  to  investigate  specificity  of  a  nonuniform  ionization  of  the  flow  we  have 
made  calculations  of  concentrations  of  charged  particles  in  plasma  sustained  by  periodically 
spaced  e-beams.  Fig.5  shows  one  of  periodical  functions  / (x)  used  in  calculations.  Space 
distribution  of  electron  concentration  produced  by  such  set  of  e-beams  is  shown  in  Fig.6.  Fig.7 
shows  another  type  of  periodical  function  /(x)  which  was  used  in  calculations.  Fig.8 

demonstrates  space  distribution  of  electron  concentration  produced  by  this  set  of  e-beams.  One 
can  see  that  space  distributions  of  electron  concentration  in  Figs. 6, 8  are  periodical  too. 

Secondly  the  nonsteady  but  uniform  approach  was  investigated  for  which  dne  +  _  jdx  =  0 


and  correspondingly  f(t,x)  =  f(t).  Results  of  calculations  presented  in  Fig.9,  corresponds  to 


0,  t  <  0 
1,  t>0 


Fig.  10  demonstrate  time-dependent  electron  concentration  in  uniform  plasma 


sustained  by  pulse-periodic  e-beam  with  a  period  T=210"5s.  And  then  the  more  complicated 
situation  in  which  concentration  of  charged  particles  varying  both  in  time  and  in  space  directions 
was  considered.  Fig.  11  demonstrates  time  and  space  dependency  of  electron  concentration  in 
plasma  sustained  by  periodically  spaced  set  of  pulse-periodic  e-beams  which  is  described  by  set 


of  equations  (1)  with  / (t,x)  =  (x) .  In  analyzing  the  results  of  calculations  the  conclusion 


was  made  that  in  case  of  both  space  and  time  dependent  plasma  concentrations  the  averaged 
concentrations  of  charged  particles  can  be  determined  by  considering  uniform  and  steady  state 
plasma  sustained  by  uniform  and  steady  e-beam  with  averaged  value  of  power  spent  on  flow 
ionization. 
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Electron  concentration  in  regions  FS,  SI  and  S2  of  MHD  controlled  inlet  was  calculated 
at  various  values  of  qum  for  flight  Mach  numbers  varying  from  6  to  12  and  various  values  of  free 
stream  dynamic  pressure.  Results  of  calculation  were  approximated  by  simple  functions  distinct 
for  different  regions  of  a  flow.  Numerical  results  and  approximation  function  for  FS  region  are 
presented  in  Fig.  12. 

Next  phase  of  this  stage  of  research  was  devoted  to  investigation  of  influence  of  induced 
electric  field  in  MHD  channel  on  electron  temperature  and  electron  concentration  sustained  by  e- 
beam  ionizer.  Results  obtained  in  Refs  [5,  6]  were  used  to  formulate  the  simple  approach  for 
calculation  of  electron  temperature  in  conditions  of  MHD  controlled  inlet.  It  is  taken  into 
account  in  this  approach  that  electron  temperature  depends  on  the  ratio  of  the  “so-called” 

E' 

effective  electric  field  E  ff=  .  to  the  gas  concentration  N,  where  E'  is  the  electric  field  in 

the  MHD  generator  in  the  reference  system  moving  together  with  a  flow,  p  is  the  Hall  parameter. 
Scheme  of  ideally  sectioned  Faraday  MHD  generator  was  considered.  It  is  shown  that  electron 
temperature  depends  on  the  magnetic  flux  density,  the  flow  velocity,  the  gas  concentration  and 
the  load  factor  k.  One  of  calculations  of  Te  in  FS  region  is  shown  in  Fig.  13.  One  can  see  that 
temperature  of  electrons  depends  significantly  upon  the  magnetic  flux  density.  Since  the  constant 
of  reactions  in  kinetic  scheme  depends  both  on  the  gas  temperature  and  electrons  temperature, 
thus  the  ionization  degree  in  the  MHD  generator  will  depend  on  the  magnetic  flux  density. 
Fig.  14  demonstrates  the  ionization  degree  in  the  MHD  generator  with  nonequilibrium 
conductivity  in  function  of  magnetic  flux  density  for  two  values  of  the  power  density  spent  on 
flow  ionization.  One  can  see  that  increasing  the  magnetic  flux  density  causes  the  ionization 
degree  to  increase.  So  effect  of  increasing  of  electrons  temperature  in  the  MHD  generator 
investigated  is  very  important  and  need  to  be  taken  into  account  in  calculations  of  a 
nonequilibrium  plasma  concentration. 

Calculated  concentrations  of  charged  particle  in  plasma  sustained  by  e-beam  in  FS,  SI 
and  S2  regions  of  the  MHD  controlled  inlet  will  be  presented  in  the  annual  report  for  various 
values  of  qion,  B,  <yx,  Mx.  Simple  formulas  for  calculations  of  important  characteristics  of  the 
nonequilibrium  plasma  will  be  presented  too. 
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x,  cm 

Fig.7.  Function  /  (x)  for  periodically  spaced  e-beams  in  form  of  a  peak  with  Dx= 4cm. 


Fig.8.  Electron  concentration  in  FS  region  for  periodically  spaced  e-beams  with  / (x)  shown  in 
fig.7  at  qion=  lW/cm3,  qm=  75  kPa,  Dx=4cm. 
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Fig.  9  Time-dependent  electron  concentration  in  uniform  plasma  sustained  by  e-beam  in  FS 

'J 

region  at  qion= lW/cm  ,  <yx=75  kPa. 

M 


Fig.  10  Time-dependent  electron  concentration  in  uniform  plasma  sustained  by  pulse-periodic  e- 
beam  with  T=2-10'5s  in  FS  region  at  t/(o;,=l  W/cm3,  75  kPa. 


Fig.  11  Time  and  space  dependent  electron  concentration  (in  cm'3)  in  plasma  sustained  by 
periodically  spaced  set  of  pulse-periodic  e-beams  with  T=110'5s  and  Dx=4cm  in  FS  region  at 
Moo=7,  <?,„„=  lW/cm\  £/x=75  kPa 
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Fig.  12  Ionization  degree  depending  on  the  parameter  qion  /N2  for  FS  region. 
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3.  Mathematical  model  to  describe  the  electron  beam  propagation  in  homogeneous 
supersonic  flow  in  presence  of  magnetic  field  will  be  developed  on  a  base  of 
Boltzmann  kinetic  equation  in  approach  of  continuous  deceleration.  Self-consistent 
approximation  with  taken  into  account  the  space  charge  and  kinetics  of  plasma 
creation  is  considered. 


One  of  the  most  important  characteristics  for  MHD  controlled  inlet  is  the  flow 
conductivity  produced  by  e-beam.  Local  concentration  of  electrons  and  of  course  local 
conductivity  of  the  flow  depends  on  the  local  value  of  the  power  density  spent  on  ionization. 
High-speed  electrons  moving  in  gas  loses  energy  in  elastic  and  inelastic  collisions  with 
molecules  of  air  it  leads  to  modification  of  electrons  energy  distribution  function.  Thus  the 
power  deposited  by  e-beam  in  the  flow  will  be  a  function  of  distance  passed  by  electrons.  So  to 
calculate  the  space  distribution  of  the  power  deposited  by  e-beam  in  the  MHD  controlled  inlet  it 
is  necessary  to  find  the  distribution  function  of  electrons  in  e-beam.  This  problem  is  solved  in  the 
stage  of  the  project.  The  second  problem  investigated  in  the  work  is  connected  to  estimation  of 
influence  of  electric  field  which  is  produced  by  uncompensated  space  charge  of  e-beam  on  the  e- 
beam  propagation. 

The  mathematical  model  for  description  of  propagation  of  electron  beam  in  homogeneous 
air  flow  in  presence  of  magnetic  field  is  developed  on  the  basis  of  Boltzmann  kinetic  equation.  In 
conditions  typical  for  MHD  controlled  inlet  the  distribution  function  for  electrons  in  e-beam 
(which  is  used  for  creation  of  nonequilibrium  conductivity),  can  be  described  by  the  Boltzmann 
kinetic  equation  for  fast  particles  moving  in  the  gas  of  slow  molecules  [1]: 


— +v— +F(iy)— =  ■ 
dt  dr  v  ’dp 


-Ncr(p)vf(r,p,Sl,t)  + 


+N  \  dSi  J  dp'v(p')-  r,  p\  Q',  t)  +  Q(  r,  p,  O,  t) 


(1) 


where  f(r,p,Sl,t)  is  the  distribution  function  of  fast  particles,  r  is  the  radius-vector  of  the  fast 
particle  position,  p  is  the  fast  particle  momentum,  fi  is  the  unit  vector  which  determines  moving 
direction  of  the  fast  particle  (Sl  =  p/p),  t  is  the  time,  v  is  the  fast  particle  velocity, 
d2a(p',£l'  — >  p,Sl) 


dpdfl 


is 


the 


differential 


cross-section  of  scattering, 


|  dpv^p ') 


d2(r(p',Sl'  ^  p,£l) 


dpd£l 


=  c(p')  is  the  total  cross-section  of  scattering,  Q(r,p,Sl,t) 


is  the  source  function  for  fast  particles,  F(r ,t)  is  the  force  acting  on  the  particle,  N  is 


concentration  of  molecules  in  gas. 

Analysis  has  shown  that  in  conditions  of  MHD  controlled  inlet  the  Larmor  radius  of 
electrons  in  e-beam  significantly  less  than  free  path  of  electrons  between  elastic  collisions.  So, 
electrons  in  e-beam  move  practically  along  the  magnetic  field.  We  will  consider  steady  state 
situation.  The  magnetic  field  is  supposed  to  be  uniform  and  perpendicular  to  the  inlet  surface. 
The  electric  field  initially  is  supposed  to  be  zero.  It  is  evident  that  in  this  case  the  energy 
distribution  function  of  electrons  will  depend  only  on  the  distance  x  from  the  surface  which  the 
e-beam  is  emerged  from.  In  this  case  the  kinetic  equation  for  fast  electrons  can  be  rewritten  in 
the  form  [2]: 

p^-  =  -Na(T)0(x,p,T)  +  iv{ dp'\ dT'a ( T')g (T\ p’  -+T,/u) <D(x, p', T)  +  Q(x, p, T )  (2) 
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where  ®(x,//,  T)  =  vf(x,ju,T)  is  the  flux  density  of  fast  electrons,  T  is  kinetic  energy  of  electron, 
ju  =  fifi# ,  fi#  is  the  unit  vector  which  determines  initial  direction  of  e-beam,  g(T',fj!  — »  T,/u) 
is  the  scattering  indicatrix  i.e.  the  probability  density  that  electron  parameters  T\fj!  take  values 
T,/u  in  result  of  scattering.  Relationship  between  //  and  //'  is  determined  by  the  formula: 

//  =  //  cos  6  +  ^/l  -  //2  sin  6  cos  xg  ,  where  9  is  the  scattering  angle  of  electron  in  the  laboratory 

system  ( cos  6  =  Qft' ),  y/-  is  the  azimuthal  scattering  angle.  In  the  general  case  cr(r')  is  the 

total  cross-section  of  electron  scattering  which  includes  both  elastic  and  inelastic  scattering.  It  is 
necessary  to  note  that  if  we  will  consider  a  nonuniform  magnetic  field  we  can  use  the  equation 
(2)  for  description  of  distribution  function  of  electrons  by  substituting  length  of  electrons  track 
along  the  magnetic  field  line  instead  of  x  coordinate. 

Kinetic  equation  (2)  is  solved  in  considering  monoenergetic  e-beam  injected  from  the 
surface  x=0  with  energy  T=To,  orthogonally  to  the  surface,  so  the  source  term  in  equation  has  a 
form: 

Q(x,  ju,  T)  =  Yb0S(x)  S(T  -  T0 )  S(ju  - 1)  (3) 

where  r/)0  is  the  initial  flux  density  in  the  electron  beam.  It  is  necessary  to  note  that  equation  (2) 
with  source  term  (3)  is  equivalent  to  the  equation  (2)  with  Q=0  and  with  boundary  condition 
9(0,M,T)  =  rtaS(T-T„)S(M-\). 

The  kinetic  equation  (2)  is  very  complicated  integro-differential  equation  which  has  not 
analytical  solution  in  the  general  case.  There  are  some  simplifications  which  allows  one  to  obtain 
approximate  solutions  of  the  kinetic  equation  [1,3].  One  of  them  is  the  continuous  deceleration 
approximation.  In  this  approximation  in  the  integral  of  collisions  instead  of  the  complex 
g(r,//->r,//)0(x,//,r)  is  used  two  first  expansion  terms  of  its  Taylor  series.  Usually  it  is 
used  in  approach  with  //  =  //'  =  1 ,  thus  function  O  is  considered  as  a  function  of  two  variables:  x 
and  T.  In  these  approaches  the  kinetic  equation  (2)  takes  the  simple  form  [1]: 

’^=-^(S(T)S>(x,T))  (4a) 

with  boundary  condition:  ®  (0,  T7)  =  rMI  ■  5 (T -T0).  Here  S(T)  is  the  so-called  stopping  power 

-  it  is  averaged  energy  lost  by  electron  with  energy  T  on  unit  length  of  its  track.  The  stopping 
power  for  electrons  in  air  can  be  obtained,  for  example,  from  Ref.  [4],  The  equation  (4a)  has  a 
simple  solution  [1]: 


dT_ 

S(T') 


(4b) 


The  space  dependency  of  the  power  density  deposited  by  e-beam  in  the  flow  in  approach 
of  continuous  deceleration  is  determined  by  relation:  ^(x)  =  ri05'(r(x))  where  f(x)  is 


solution  of  equation:  |  dT'/S  (T'j  =  x .  The  approach  considered  doesn’t  take  into  account 

T(x) 

fluctuation  of  energy  losses  and  dispersion  of  electrons  ranges  [1]  and  of  course  it  can  be  used 
only  for  simple  estimations. 

In  order  to  obtain  analytical  solution  of  equation  (2)  more  realistic  than  solution  (4b)  it  is 
necessary  to  have  simple  expression  for  the  scattering  indicatrix  g  (T', //'  — that  really 


describes  modification  of  electron  energy  in  collisions  with  molecules  of  a  flow  in  a  wide  range 
of  electrons  energies.  We  have  analyzed  differential  cross-sections  which  are  used  in  modeling 
elastic  and  inelastic  scattering  of  fast  electrons  in  gases  and  solids.  The  main  channel  of  electron 
energy  degradation  is  inelastic  scattering.  When  the  energy  loss  of  fast  electron  in  inelastic 
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collisions  with  atom  of  a  medium  is  significantly  higher  then  the  binding  energy  of  electrons  in 
the  atom,  the  inelastic  collisions  can  be  approximately  considered  as  collisions  with  free 
electrons.  In  this  case  there  are  simple  equations  for  description  of  the  differential  cross-section, 
see  e.g.  [2]: 

dcr  _  27rr^mec2  1 
dW  ~  p2  W2 

where  W  is  the  energy  lost  by  the  fast  electron  in  collision  with  a  motionless  electron,  ro  is  the 
classical  electron  radius,  me  is  electron  mass,  c  is  velocity  of  light,  t  =  T  /mec2 ,  P  =  v/c . 

When  W  is  compared  with  the  binding  energy  it  is  necessary  to  use  more  complicated 
expressions  for  the  differential  cross-section,  see  e.g.  [2], 

In  the  project  we  have  developed  several  modeling  differential  cross-sections  for 
inelastic  scattering  of  fast  electrons  on  air  molecules  in  form: 


1  +  - 


W2 


-  +  - 


w_ 

\2  rri  2 


2r  +  l 


W 


(: T-W)  (r  +  l)  T  (r  +  l)z  T-W 


(5) 


dcr 

~dW 


27irpnec 


P 


■0(T,W) 


(6) 


Evident  requirement  for  function  &(T,W)  is  providing  good  agreement  with  known 
values  of  the  stopping  power.  In  our  designations  the  stopping  power  [4]  is  described  by 
equation: 


IV„ 


S(r)  =  NZ-^3£  |  W-0(T,W)dW 

P  w, . 


(7) 


Here  we  consider  several  functions  d>(7’,IE)  which  were  proposed  in  our  work: 


1  =  Wmin=tn(T),  Wmax=T/2. 

W 

2.  < D2(r,  =  Wmin=tl2(T),  Wmax=T. 

3.  03 <T,  W)  =  - - - - - ;  Wmi„=0,  Wmax=T/2. 

(^  +  *73^)) 

4.  O 4(T,  W )  =  - - - ^ ;  Wmin=0,  Wmax=T. 

(W  +  ilA{T)) 


These  approximations  allow  one  to  obtain  simple  analytical  expressions  for  stopping  power.  In 
using,  for  example,  QX(T,W)  we  obtain  j  W ■Ol(T,W)dW  =  In Functions  ^(T), 


jj2  (r) ,  //3  (T)  and  t]4  (T)  are  fitting  functions,  which  are  chosen  as  solution  of  equation  (7)  in 

using  values  of  S(T )  from  Ref.  [4],  Fig.l  shows  dependencies  of  ( T )  and  /;3  ( T )  from  energy 

of  fast  electron.  One  can  see  that  these  fitting  functions  are  monotonic  and  can  be  approximately 
described  by  simple  power  functions. 

Calculations  have  shown  that  a  value 


x=zj 


dcr  2nr2m  c 2  W° 

dW  =  Z-  0  e 


tVmi 


dW 


P 1 


j  0(T,W)dfV  =  Z 


2/rr0mec 


Wm, 


p2 


t(t), 


is  slow  depending  function  of  electron  energy.  Taking  into  account  this  result  we  propose 
additional  way  to  choose  fitting  function  from  condition  £  =  const  -  £ .  In  such  approach  we 
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have  obtained  e.g.  ;/,  =  7/(2  +  4.87  -107  r/?2  j .  The  value  X  can  be  considered  as  a  total  cross- 

section  of  inelastic  scattering  of  fast  electron  on  molecule.  Fig.2  shows  the  stopping  power 
calculated  under  equation  (7)  in  using  ®  =  ®t .  Three  approaches  were  used  in  calculations  of 

rjx  (7) .  One  can  see  that  all  approaches  are  in  good  accordance  with  the  set  of  data  from  [4], 

Distinguish  of  analytical  dependencies  in  the  range  of  low  energies  7<0.3keV  is  not  significant 
in  investigations  of  space  distribution  of  power  deposited  by  high  energy  electrons. 

In  order  to  make  next  step  to  solution  of  kinetic  equation  we  have  to  determine 
probability  density  g(T'  — >  7)  that  electron  changes  its  energy  from  7'  value  to  7  in  result  of 

scattering  on  molecule.  In  particular  in  using  ®  =  ®4  we  obtain: 


g(T'^T) 


^4(7')-(7'  +  ^4(7')) 
7'-(7'-7  +  t/4(7'))2 


(8) 


The  angle  of  scattering  6  can  be  calculated  in  using  the  conservation  laws.  Estimations 
show  that  the  angle  of  scattering  of  electrons  at  typical  for  condition  of  MHD  controlled  inlet 
energies  is  very  small.  Thus  in  this  stage  of  project  we  restrict  our  investigations  by  considering 
situation  with  cos  0  =  1  and  //  =  //.  In  these  approaches  we  have  obtained  solution  of  kinetic 
equation  (2)  in  form  of  expansion  on  the  number  of  collisions: 

oo 

®(2,^,7’)=rM- £jr,(z)F,(r)5(/,-i)  (9) 

k= 0 

where  z  =  NTx ,  functions  Xk  (z)  and  Fk{T)  satisfies  to  recurrence  relations: 


dz 


+Xl(z)  =  Xi_,(z) 


T 

F,(r)=jg(r->7')-F1.,(r)rfr 


(10) 


T 

with  X0  (z)  =  exp(-z)  and  F{)(T)  =  S(T -T{)) .  Equation  for  Xk  (z)  can  be  explicitly  solved: 
Xk  (z)  =  zk  exp {-z)jk\  at  k  >  0 .  The  obtained  solution  allows  one  to  calculate  energy  spectrum 
for  electrons  at  any  given  point  of  z  or  x  =  z/ NT.  by  the  equation  F  (T)  =  ^Fk(T)  Xk  (z)  .  Also 

k 

we  can  calculate  averaged  energy  of  electrons  depending  on  the  z  coordinate: 


T  (4  ■ =  2X  Wjr-F,  (T)dT  =2X  (z)  Tt 

k  o  £ 


(11) 


where  Tk  is  averaged  energy  of  electrons  after  k  collisions.  And  of  course  we  can  calculate  space 

distribution  of  the  energy  lost  by  electron  in  unit  of  length  -dT/dz  (or  - dT/dx )  value  by  the 
relation: 


dT  v  ~  dX^  (z)  y,  - 
dz  4"  *  dz  k  k 


[zk  -kzk  1  j 


k\ 


exp(-z) 


(12) 


Space  distribution  of  the  power  density  deposited  by  e-beam  in  flow  can  be  determined  in 
these  designations  by  the  relation  q  =  r/)(l  •  (- dT/dx ) . 


Fig.3  shows  comparison  of  energy  spectrums  for  electrons  calculated  in  using  two 
approaches:  ®  =  ®2  and  ®  =  ®4.  One  can  see  that  the  spectrums  practically  coincide  despite 
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distinguish  of  functions  g(T'  — >  T)  in  these  approaches.  This  coincidence  confirms  the  criterion 

of  choice  for  the  modeling  differential  cross-section  used  in  our  work.  Fig.4  demonstrates 
comparison  of  the  energy  spectrum  for  e-beam  transmitted  through  carbon  with  a  thickness  of 
230A  which  was  computed  by  Monte-Carlo  method  in  [5]  and  calculated  in  our  work  in  using 
approach  with  O  =  d>4 .  Fitting  function  ;/4  (T)  calculated  for  carbon  target  was  used  in  this  case. 

A  good  accordance  of  these  results  demonstrates  validity  of  the  method  of  simplification  for 
description  of  collision  of  fast  electrons  with  atoms  and  molecules  suggested  in  our  work.  Fig.5 
demonstrates  the  energy  spectrum  in  e-beam  with  initial  energy  T0=10keV  at  various  values  ofx. 
One  can  see  that  increase  of  the  length  the  e-beam  pass  through  decreases  the  maximal  energy  in 
e-beam  and  widens  the  energy  spectrum. 

Fig.6  demonstrates  the  space  distribution  of  the  energy  lost  by  electrons  in  unit  of  length. 
Calculations  in  two  approaches  are  compared  here.  One  can  see  that  continuous  deceleration 
approach  provides  agreement  with  the  new  developed  approach  (Equation  (12))  only  at  small 
values  of  coordinate  x.  Fig.7  demonstrates  dependencies  of  dT/dx  for  e-beams  with  various  initial 
energies  To.  One  can  see  that  increase  of  To  leads  to  increasing  the  length  of  e-beam  penetration 
into  the  air  flow. 

In  order  to  investigate  influence  of  electric  field,  which  is  produced  by  uncompensated 
space  charge  of  e-beam,  on  the  e-beam  propagation  we  need  to  know  dependency  of  the  e-beam 
flux  density  on  the  spatial  coordinate  -  Tb  (x)  or  T,,  (z).  In  using  solution  of  Boltzman  kinetic 

equation  in  form  (9,  10)  we  can  introduce  rft  (z)  by  the  relation: 

1  To  oo  ^0 

rs  (2)  =  fd/t J  9(z,fi,T)dT  =  rsa  -2X  (z).  J  Fk  (T)dT  (13) 

-1  T,  k= 0  t, 

1  1th  1th 


where  Tth  is  threshold  energy.  This  formulation  implies  that  in  case  when  energy  of  electron  is 
less  than  threshold  energy  ( T  <  T,h  X  it  is  removed  from  e-beam.  This  process  can  be  interpreted 
as  e-beam  thermalization.  Fig. 8  demonstrate  dependency  of  rfi (x)  at  Tth=  leV.  System  of 

equations  used  in  our  analysis  for  description  of  forming  the  electric  field  E(x)  in  the  track  of  e- 
beam  propagation  is  quite  similar  to  one  used  in  Ref.  [6]. 


dfle  ,  T  r3  dYb 

+  ion  ~  PeM  ~  ~T 

ot  ox  ax 


PeM 


ion  r'  ei  e  i 


dn,  ,  dr,  _7 

d  t  dx 

dE  e  /  \ 


(14) 


(3/2  0/2 

wherere(x)  =  -//e«e(x)£,(x)-Z)e — -,  T;  (x)  =  //  h  (x)  £ (x)  -  Dt — L ,  //,  and  //,.  are  electron 

dx  dx 

and  ion  mobility,  De  and  D  are  electron  and  ion  diffusion  coefficients,  Iion  is  the  ionization  rate, 
ne  is  the  electron  concentration,  n,  is  the  ion  concentration,  nb  is  concentration  of  fast  electrons  in 
e-beam  ( nb  =  f/;  (x)/v  (x) ),  /T  is  the  electron-ion  dissociative  recombination  rate  coefficient. 


Ionization  rate  is  determined  by  the  relation  Iion  = 


rM> -(-dT/dx)  _  jb -(-dT/dx) 


W,„ 


eW;„ 


,  where  jb  is  the 


current  density  of  e-beam,  Wion  is  ionization  cost.  In  equations  (14)  we  consider  simplified 
kinetic  of  ionization  (negative  ions  are  not  included  in  this  kinetic  scheme)  because  the  main  task 
of  this  stage  of  analysis  is  estimation  of  role  of  the  electric  field  E(x)  on  propagation  of  the  e- 
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beam.  If  conditions  at  which  this  electric  field  significantly  modifies  e-beam  passage  in  the  air 
flow  will  be  found  the  more  complicated  scheme  for  corresponding  calculations  will  be  realized. 
We  consider  steady  state  solution  for  which  dne /dt  =  dnjdt  =  0 .  Values  of  electron  mobility  and 
diffusion  coefficient  were  calculated  according  to  data  from  [7]  in  supposition  that  there  is 
transversal  electric  field  in  MHD  channel  of  the  MHD  controlled  inlet  £’1=10V/cm.  To  estimate 
influence  of  electric  field  E(x)  on  the  e-beam  moving  in  the  air  plasma  we  have  calculated  a 

Xmax 

value  of  potential  drop  AU  in  the  region  the  e-beam  passing  through:  A U  =  -  J  E  (x)dx . 

0 

Relative  value  8  =  \eAU\/T0  determines  relative  effect  of  electric  field  E(x)  on  additional 

deceleration  of  e-beam.  Calculations  for  typical  conditions  for  MHD  controlled  inlet  have  shown 
that  parameter  8  <  0.01 .  So  effect  of  additional  deceleration  of  e-beam  by  this  electric  field  can 
be  neglected. 

In  the  next  stage  of  the  project  the  analytical  solution  of  Boltzmann  kinetic  equation  will 
be  generalized  to  the  case  of  inhomogeneous  air  flow  in  which  the  flow  density  will  be  a 
function  of  x  coordinate.  Effect  of  changing  the  moving  direction  of  electrons  in  scattering  will 
be  accounted  in  calculation  the  spatial  distribution  of  the  power  deposited  by  e-beam.  In 
calculations  of  potential  drop  AU  in  the  region  the  e-beam  passing  through  the  self-consistent 
approach  will  be  considered  taking  into  account  that  electric  field  E±  induced  in  the  channel  of 
the  MHD  generator  depends  on  the  flow  and  the  magnetic  field  parameters. 


CONCLUSIONS 


In  the  third  stage  of  the  project  modeling  differential  cross-section  for  inelastic  scattering 
of  fast  electrons  on  atoms  and  molecules,  which  ensure  good  describing  the  stopping  powers, 
were  offered.  In  using  the  suggested  differential  cross-sections  we  have  obtained  new  analytical 
solution  of  Boltzmann  kinetic  equation  for  distribution  function  of  fast  electrons  moving  in 
homogeneous  air  flow.  New  solution  was  compared  with  approach  describing  e-beam  passing  in 
gases  in  approximation  of  continuous  deceleration.  It  is  shown  that  continuous  deceleration  may 
be  used  to  calculate  the  power  deposition  by  e-beam  only  in  initial  part  of  its  track.  Analysis  of 
influence  of  electric  field,  which  is  produced  by  uncompensated  space  charge  of  e-beam,  on  the 
e-beam  propagation  was  made. 
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Fig.  1  Fitting  functions  (T )  and  r]}  (T ) 


T,  keV 

Fig.2  Stopping  power  for  electrons  in  air  calculated  in  using  O  =  O,  at  various  choices  of  fitting 
function  77,  (T). 
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Fig.3.  Energy  spectrum  for  an  e-beam  in  air  at  x=2cm.  To=2keV,  N=1017  cm"3 


T,  eV 

Fig.4.  Energy  spectrum  for  an  electron  beam  transmitted  through  carbon  with  a  thickness  of 
230A  with  T0=1600  eV. 
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T,  keV 

Fig.5.  Energy  spectrum  for  an  e-beam  in  air  at  various  values  of  x.  To=2keV,  N=1017  cm'3 


x,  cm 

Fig.6.  Space  distribution  of  the  energy  lost  by  electrons  in  unit  of  length  calculated  in  two 
approaches. 
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x,cm 

Fig.7.  Space  distribution  of  the  energy  lost  by  electrons  in  unit  of  length  for  two  values  of  initial 
energy  of  e-beam,  calculated  under  equation  (12). 


x,  cm 

Fig.8.  Spatial  dependency  of  e-beam  flux  density  at  Eth=leV 
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4.  Mathematical  model  to  describe  the  electron  beam  propagation  in  nonuniform 
supersonic  flow  in  presence  of  magnetic  and  electric  fields  will  be  developed  on  a 
base  of  Boltzmann  kinetic  equation  in  approach  of  continuous  deceleration.  Self- 
consistent  approximation  with  taken  into  account  the  space  charge  and  kinetics  of 
plasma  creation  is  considered.  Space  distribution  of  electron  and  ion  concentrations 
will  be  calculated  in  typical  for  MHD-controlled  inlet  conditions 


In  the  stage  of  the  project  we  have  investigated  e-beam  propagation  in  general  case  of 
nonuniform  supersonic  flow.  Influence  of  space  charge  produced  by  e-beam  in  plasma  on 
propagation  of  e-beam  and  on  formation  of  spatial  distribution  of  electron  and  ion  concentrations 
in  the  plasma  is  analyzed  here  too. 

At  first  we  consider  main  results  obtained  in  this  stage  in  analysis  of  space  charge  and 
electric  field  distribution  in  track  of  e-beam  propagation.  System  of  equations  used  in  our 
analysis  coincides  with  the  system  considered  in  previous  stage  of  the  project: 


dn„  8F, 


-  +  - 


dt  dx 


-  -  /• 
1ion 


■fie 


dx 


dm  dV: 

— L  +  — L  =  7, 
dt  dx 


Peineni 


(1) 


dE_ 

dx 


-{ni-ne-nb) 

^0 


Main  designations  used  here  are  conventional.  Coordinate  x  is  directed  along  the 
propagation  direction  of  e-beam  and  measured  from  the  surface  the  e-beam  is  emerged  from,  Im 
is  the  ionization  rate,  rib  is  concentration  of  fast  electrons  in  e-beam. 

We  take  into  account  that  electron  mobility  and  the  dissociative  recombination  rate 
coefficient  Pei  are  functions  of  electron  temperature  which  according  to  results  obtained  in  the 
second  stage  of  the  project  depends  on  the  effective  electric  field  Eeff  in  channel  of  MHD 
generator.  The  electric  field  in  this  case  will  be  a  combination  of  a  longitudinal  electric  field 
(oriented  along  the  magnetic  field  lines)  produced  by  uncompensated  space  charge  of  e-beam, 
and  the  transverse  electric  field  E±  induced  in  the  MHD  generator.  In  the  stage  of  the  project  we 
have  obtained  analytical  solution  of  equations  (1)  in  justified  neglecting  both  the  ion  flux 
Fj  «Teand  the  gradient  of  the  ion  flux  dTf/dx  « Iion.  In  assuming  rib  =  0  the  solution  of 


equations  (1)  for  electric  field  takes  on  the  form: 


E(x)  = 


X 


max 


x' 


2  |  5(x')exp(-|24(x")<7x")d'x' 

X  X 


(2) 


where  A(x)  =  ^lionPpPj^e  =  e^(x)  ;Xmax  js  maximal  length  of  penetration  of  e- 
£oPeirb(x) 

beam  into  a  flow  (the  so-called  range  of  incident  electrons).  Considering  the  real  values  of  rib(x) 
insignificantly  increase  the  electric  field.  The  relative  increment  of  E(x)  is  typically  less  than 
0.1%  and  it  can  be  neglected. 

The  solution  (2)  allows  us  to  estimate  the  influence  of  the  electric  field  E(x)  on  the 
penetration  of  e-beam  in  the  air  plasma  in  a  wide  spectrum  of  conditions,  taking  into  account 
self-consistency  of  the  problem.  A  value  of  potential  drop  .  U  in  the  region  the  e-beam  passing 


through  A U  =  -  |  E(x)dx  determines  what  energy  the 
0 


e-beam  loses  not  in  collisions  but  in 
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deceleration  in  the  electric  field.  Relative  value  S  =  \eAU\/To  determines  relative  effect  of 

electric  field  E(x)  on  additional  deceleration  of  e-beam.  In  consideration  of  results,  obtained  in 
previous  stage  of  the  project  (dependencies  for  Iion(x)  and  T^(x)),  we  have  obtained  simple 
relations  which  allow  us  to  estimate  the  value  8 : 


8<8  =  ■ 


l ejbPeiWion  /-hmax  \3/2 


(3a) 

where  jb  is  the  current  density  in  the  e-beam,  Wion  is  the  ionization  cost.  In  using  results  of 
[1],  which  allow  us  to  determine  experimentally  measured  values  of  xmax  in  terms  of  To  and  N, 
the  relation  (3  a)  can  be  rewritten  in  the  form: 


^  =  0.74.10-4-^SZ|^V  • t 06 

z  t  AT  J 


(3b) 


where  (3  = - — - — ,ju  = - — .  Parameter  L  is  required  length  of  the  flow 

2-10 cm5  /  s  io5  cm1/ 

/V  • s 

region  in  which  the  e-beam  produces  conductivity.  Dimensions  of  parameters  jb,  N  and  L  in 
relation  (3b)  are  correspondingly  mA/cm2,  cm'3  and  cm.  Fig.l  shows  that  increasing  the  Eeff/N 

value  leads  to  increasing  the  8  parameter.  Increasing  the  air  density  N  causes  the  parameter  8  to 
be  increased  insignificantly.  It  is  evident  from  Fig.l  that  in  conditions  typical  for  MHD 

controlled  inlet  the  parameter  8  « 1 .  Thus  in  these  conditions  the  influence  of  electric  field  on 
propagation  of  ebeam  in  the  flow  is  insignificant  and  it  can  be  neglected.  So  our  investigations 
show  that  the  spatial  distribution  of  the  energy  deposited  by  e-beam  in  conditions  of  MHD 
controlled  inlet  practically  doesn’t  depend  on  the  electric  field  generated  in  the  region  of  the  e- 
beam  propagation.  Spatial  distribution  of  concentrations  of  charged  particles  and  conductivity  of 
flow  sustained  by  e-beam  will  be  functions  of  flow  parameters  and  e-beam  characteristics. 

In  the  previous  stage  of  the  project  it  was  shown  that  the  flux  density  ®(x,//,  T)  of  fast 
electrons  in  e-beam,  which  allows  one  to  determine  the  spatial  profile  of  the  energy  deposited  by 
the  e-beam  in  a  flow,  is  described  by  the  Boltzmann  kinetic  equation  in  the  form: 

V^  =  -Na(T)®(x,Ju,T)  +  NfdJu'fdT'(j(r)g(r,Ju'^T,JuYt>(x,M',T')  +  Q(x,Ju,T)  (4) 

Solution  of  the  equation  was  found  in  considering  g(T',/u'— >  T,/u)  =  g(T'^>  T )  •  8(ju'-ju ) 
and  in  assuming  that  a  flow  concentration  N  doesn’t  depend  on  the  x  coordinate.  Here  8{x)  is 
Dirac  delta-function.  In  conditions  of  MHD  controlled  inlet  a  supersonic  flow  usually  is 
nonuniform  and  a  concentration  N  depends  on  the  x  coordinate,  i.e.  N=N(x).  Here  we  generalize 
the  solution  obtained  in  previous  stage  of  the  project  to  situation  with  N=N(x),  and  develop  a 
recipe  of  calculation  of  spatial  profile  of  the  energy  deposition  by  e-beam  in  this  case.  Let  us 
introduce  a  new  variable  x  by  the  relation: 


X 


0 

where  No  is  characteristic  concentration.  In  using  this  variable  the  equation  (4)  takes  on 
the  form: 


Ju—  =  -N0cr(T)®(x,jU,T)  +  N0jdJu'jdrcT(r)g(T',ju'^>T,ju)®(x,Ju',r)  +  Q(x,jU,T)  (6) 

One  can  see  that  equation  (6)  coincides  with  the  equation  for  uniform  flow  in  case  of 
N=No.  So  solution  of  kinetic  equation  (4)  in  general  case  of  N=N(x)  can  be  obtained  from 
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solution  of  equation  (6)  by  using  evident  transformation:  <£»(*,  u,T)  -  d>(x  =  f  dx,u,T). 

0  N” 

Spatial  distribution  of  energy  deposition  dT/dx  in  the  nonuniform  flow  can  be  expressed  through 
dT/dx  by  the  relation: 


dT  dT  vx  dx  dT 

—  =  — (x  =  x(x)) - =  — 

dx  dx  dx  dx 


( 


=  1 


N(x) 

Nn 


dx 


N(x) 

iVn 


(7) 


V  o  “  j 

Here  the  function  dT/dx  is  calculating  through  ®(x,//, T)  in  accordance  with  results 
obtained  in  previous  stage  of  the  project.  It  is  easy  to  see  from  equation  (7)  that  dependency  of  a 
flow  concentration  on  the  x  coordinate  significantly  influences  on  the  profile  of  the  energy 
deposited  by  e-beam  in  the  flow.  In  particular,  for  step  function  of  N(x)  caused  by  oblique  shocks 
in  a  supersonic  flow  the  dT/dx  profile  will  be  a  step  function  too.  The  obtained  results  allow  one 
to  calculate  the  spatial  profile  of  the  energy  deposited  by  e-beam  in  a  nonuniform  flow  when  the 
spatial  profile  of  gas  concentration  in  the  flow  is  known.  Figs.  2-5  demonstrate  that 
nonuniformity  of  a  flow  significantly  influences  on  the  profile  of  the  energy  deposited  by  e-beam 
in  the  flow.  Results  obtained  here  and  in  previous  stages  of  the  project  give  us  opportunity  to 
calculate  spatial  distributions  of  concentration  of  charged  particles  sustained  by  e-beam  in  a 
nonuniform  flow  of  MHD  controlled  inlet.  Figs. 6, 7  demonstrate  distributions  of  concentrations 
of  electrons  and  negative  ions  which  were  calculated  for  nonuniform  flow  modeled  by  means  of 
step  function  Ni(x),  shown  in  Fig.2.  One  can  see  that  concentrations  of  charged  particles 
significantly  change  along  the  e-beam  track. 

The  next  step,  made  in  this  stage  of  the  project,  consists  in  obtaining  more  precise 
solution  of  the  Boltzmann  kinetic  equation  (4)  which  takes  into  account  that  in  inelastic 
scattering  not  only  electron  energy  is  changed  but  also  deflection  of  electron’s  path  occur. 
According  to  [2]  the  cosine  of  the  angle  of  scattering  of  fast  electron  in  inelastic  collisions  with 

molecules  can  be  calculated  under  the  formula:  cos©  =  *J(1-w)(t  +  2)/((1-w)t  +  2)  ,  where 
w=  W/T,  W  is  the  energy  lost  by  the  fast  electron  in  collision  with  an  electron  in  molecule  of  gas 

9 

(W  =  T-  T’),  T’  is  the  energy  of  electron  after  collision,  r  =  T  /  mec  ,  me  is  electron  mass,  c  is 

velocity  of  light.  Average  value  of  cos  0  for  an  electron  with  energy  T  is  determined  by  the 
relation: 


l 

cos #(/')=  jcos 6  g(T  ->T')dT' 


(8) 


Fig.8  shows  dependency  of  average  value  of  cosine  of  the  angle  of  scattering  of  fast 
electron  in  inelastic  collisions  upon  the  energy  of  electron.  In  calculations  we  use  function 
g(T  — >7”)  proposed  in  previous  stage  of  the  project.  One  can  see  that  at  T>lkeV  the 

cos  ©magnitude  is  very  close  to  1.  But  in  the  project  we  show  that  a  moving  direction  of  fast 
electron  significantly  changes  in  result  of  multiple  collisions  with  molecules.  It  is  demonstrated 

in  Fig.9,  where  cos  ©is  average  value  of  cosine  of  the  angle  of  scattering  of  a  fast  electron  in 
result  of  an  inelastic  collision  number  k,  and  ///f  is  average  value  of  cosine  of  angle  between  the 

initial  moving  direction  of  e-beam  and  moving  direction  of  electrons  of  the  e-beam  changed  in 
result  of 

k  collisions.  Parameter  To  is  initial  energy  of  electrons  in  the  e-beam.  Value  of  cos  *0  is 
determined  bv  the  eauation: 

cos  6k  =  Jcos  1)  •  FkA(T)dT  ^ 

0  _ 

where  Fk  (7  )  is  the  energy  distribution  function  of  electrc  //.  in  e-beam  after  k  collisions, 
which  was  calculated  in  previous  stage  of  the  project.  Transformation  of  value  in  result  of 
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k 


multiple  inelastic  collisions  is  determined  by  the  equation:  jU/(  =  jUq  ■  ]^[  cos  ©y  .  In  our 

7=1 

investigations  we  assume  po=l  which  corresponds  to  the  e-beam  directed  along  the  x  axis.  It  is 
evident  that  value  of  « 1  corresponds  to  electrons  which  are  moving  practically  without 

deviations  from  the  initial  trajectory.  Value  of  «  0  corresponds  to  practically  isotropic 
movement  of  electrons.  It  is  important  to  note  that,  according  to  results  shown  in  Fig.9,  at 

n  -  - 

k  «  10  we  have  cos©/f  «  1  and  «  0 .  So  there  are  conditions  at  which  the  moving  direction 
of  electrons  practically  does  not  change  in  single  scattering  but  their  movement  after  k  collisions 

becomes  practically  isotropic.  Fig.10  demonstrates  that  at  k « 10  the  average  energy  of 

T 
±  0 

electrons  7^  =  Jr  ■  F/<(T)dT  is  large  enough.  So  we  show  that  in  order  to  correctly  calculate  the 
0 

spatial  profile  of  the  energy  deposited  by  e-beam,  it  is  necessary  to  take  into  account  deflection 
of  electron’s  path  in  collisions. 

In  order  to  obtain  analytical  solution  of  equation  (4)  we  assume  that  the  scattering 
indicatrix  g(7”,//'-»  T,ju)  has  a  form:  g(T',/u'^T,ju)  =  g(T'-*T)-g(ju'^>  .  This 

supposition  implies  that  change  of  electron  energy  in  a  single  collision  is  independent  of  the 
change  of  a  moving  direction  of  electron  in  this  collision.  This  approach  is  absolutely  correct  at 
the  initial  part  of  electrons  track  where  Ju^  « 1  and  at  the  ending  part  of  the  track  where  /7/c  «  0 
[3],  As  a  function  g(T'->  T)  we  use  one  of  functions  developed  in  previous  stage  of  the  project 


w„. 


in  approach  for  which  S  =  Z  j  dW  =  S  is  constant  value. 


w 


dW 


In  this  case  kinetic  equation  (4)  can  be  rewritten  in  the  form: 

M™=  -«*(*,  M,  T)+  J  dfi’gifi’  -+  //)f  </n?(r  ->  T)a>(z,ju’,T')+ 

CZ  -I  7  ? 

■tVlW-W"1) 


(10) 


where  dimensionless  coordinate  z  =  J  N(x'  )I.dx'  take  into  account  dependency  of  a  flow 

0 

concentration  upon  x  coordinate.  The  source  term  Q(z,  ju,  T )  =  Y ^6(2)3 (T  -  Tq)S(/u  - 1) 

in  equation  (10)  implies  that  an  e-beam  with  initial  energy  T=Tf,  is  emerged  from  the  surface  z=0 
along  the  x  axis  (//  =  1) .  It  is  easy  to  show  that  solution  of  equation  (10)  has  a  form: 

OO 

®(z,Ju,T)  =  rb0YJ  V/C(z,//)  •  F/^T).  By  substituting  this  expansion  in  equation  (10)  we 
k= 0 

obtain  the  set  of  recurring  equations  for  functions  X/C(z, //)  andF^(T)  : 

dXk{ztfi) 


M- 


■  +  Xk  (z,n)  -  J  g(fi'->u)Xk  ,  (z,fi’)dn’ 


Fk{T)=\g(r^T)Fk_x{r)dr 


(ii) 


X0  (z,  fi)  -  exp(-z)  ■  (//  - 1) ,  xk(0,ft)  l>0=0r  Xk{z,M)  Ut^O,  F0{T)  =  S(T -T0) . 

Equations  for  energy  distribution  function  Fk(T)  are  solved  in  previous  stage  of  the 
project.  In  this  stage  of  the  project  we  obtain  solution  of  equations  (11)  for  functions  X/(  (z,  f.i ). 
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To  solve  equations  (11)  analytically  we  propose  simple  approximation  for  g( //'  — >  //)  which  is 
similar  to  the  so-called  transport  approach  [4]: 

g(M'  ->  /')  =  gk  (m'  -»//)  =  cos  9k  ■  5 (//  -  ju)  + 1  ■[£(//-//)  +  £(//  +  //)]  (12) 

Here  the  term  +  c>(//  +  //)]/ 2  models  the  isotropic  scattering  in  using  value  of 

ju"=l/yf3  .  The  indicatrix  (12)  provides  correct  description  of  features  of  scattering  of  electrons 
depending  on  the  number  of  collisions  k.  According  to  Fig.9  and  equation  (12),  electrons 
practically  do  not  change  the  moving  direction  at  small  values  of  k,  but  in  increasing  value  of  k 
the  scattering  of  electrons  becomes  more  isotropic.  In  the  project  it  is  shown  that  the  indicatrix 
(12)  guarantees  accurate  description  of  dependency  of  Ji \  upon  k. 

In  the  stage  of  the  project  we  have  obtained  solution  of  equations  (1 1)  in  the  form: 

Xk(z^l)  =  Ak{z)-5{^-\)+Bk{z)-5{i,i-f;i)  +  Ck{z)-5{M  +  ^  (13) 

where 

k  _ _  k  _  k 

4  (-)  =  FI cos 4  •  exp (-2)  ■  Bk  ~ exp (-z ) 

B,  ( z)  =  exp  (---)y  bfj  ■  zJ  +  exp  (-  z/ ft)  V  bk2  j  ■  z1 

J= 0  ;=0 

Q  (-)  =  exp(-z)  •  zy  +exp(-z//7)y4J  ■ 

/=0  j=0 

Expressions  for  coefficients  bfj  and  cfj  are  bulky  enough  and  will  be  presented  in  the 

final  report.  The  obtained  solution  allows  us  to  calculate  the  flux  density  function  <f >(z,//,T)  and 

average  value  of  energy  of  e-beam.  In  using  this  solution  the  average  energy  of  electrons  in  e- 

00 

beam  is  determined  by  the  equation:  T(z)=  ^  y  •  \A/C  (z)  +  jj  -(B/c  (z)  -  C/c  (z))] .  The  spatial 

k= 0 

distribution  of  the  energy  deposited  by  e-beam  in  a  flow  is  determined  by  the  equation: 


d T 
dz 


(H) 


Fig.  11  demonstrates  the  spatial  distribution  of  the  energy  deposited  by  e-beam  in  air 
calculated  under  equation  (14)  (red  curve)  in  comparing  with  one  calculated  in  previous  stage  of 
the  project  (black  curve).  One  can  see  that  the  dependencies  differ  significantly.  It  demonstrates 
that  the  effect  of  deflection  of  electrons  path  in  inelastic  collisions  is  very  important  and  it  needs 
to  be  taken  into  account  in  calculations  of  the  energy  deposited  by  e-beam  in  a  flow. 

In  the  next  stage  of  the  project  the  main  processes  which  are  responsible  for  energy  and 
spatial  transformation  of  e-beam  propagating  in  the  MHD  controlled  inlet  will  be  analyzed.  The 
algorithm  for  Monte  Carlo  simulation  of  e-beam  propagation  in  nonuniform  supersonic  flow  will 
be  developed.  Monte  Carlo  code  for  calculation  of  e-beam  propagation  in  the  MHD  controlled 
inlet  will  be  developed.  Testing  the  Monte  Carlo  code  will  be  done. 


CONCLUSIONS 


In  the  fourth  stage  of  the  project  the  analytical  solution  of  the  Boltzmann  kinetic  equation 
is  generalized  to  the  case  of  nonuniform  supersonic  flow  in  which  the  flow  density  is  a  function 
of  x  coordinate.  In  calculations  of  potential  drop  A  U  in  the  region  of  the  e-beam  passage  the 
self-consistent  approach  is  considered.  It  is  shown  that  the  spatial  distribution  of  the  energy 
deposited  by  e-beam  in  conditions  of  MHD  controlled  inlet  practically  doesn’t  depend  on  the 
electric  field  produced  by  uncompensated  space  charge  in  the  region  of  the  e-beam  propagation. 
Spatial  distributions  of  concentrations  of  charged  particles  sustained  by  e-beam  in  a  supersonic 
flow  are  calculated  for  typical  conditions  for  MHD  controlled  inlet.  Effect  of  changing  the 
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moving  direction  of  electrons  in  scattering  is  accounted  in  calculation  the  spatial  distribution  of 
the  power  deposited  by  e-beam.  It  is  shown  that  the  effect  of  deflection  of  electrons  path  in 
inelastic  collisions  is  very  important  and  it  needs  to  be  taken  into  account  in  calculations  of  the 
energy  deposited  by  e-beam  in  a  flow. 


Fig.  1 .  Parameter  Jin  MHD  controlled  inlet  at  various  concentrations  of  air,  shown  near  the 
curves:  jb=lmA/cm2,  L=100c  m. 
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Fig. 2  Distributions  of  concentration  N(x)  used  in  calculations  of  the  energy  deposition  profiles. 


Fig.  3,  Spatial  distribution  of  energy  deposited  in  air  by  e-beam  with  To=10keV,  No=1017ctn  1 
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Fig.4  Distributions  of  concentration  N(x)  used  in  calculations  of  the  energy  deposition  profiles. 


Fig,  5,  Spatial  distribution  of  energy  deposited  in  air  by  e-beam  with  T0=10keV,  No=1017cm 
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Fig,6,  Concentration  of  electrons  in  nonuniform  flow  with  N  Nj(x)  (see  Figr2)  at  No=1017cm  3, 
sustained  by  e-beam  with  T0=10keV  and  jb=lm  A/cm2,  EegNo^51d, 
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Fig. 7.  Concentration  of  negative  ions  in  nonuniform  flow  with  N=Nj(x)  (see  Fig.2)  at 
No=10l7cm'3,  sustained  by  e-beam  with  To=10keV  and  jh=lmA/cm2,  EetfNo^STd. 
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Fig.  8.  Average  cosine  of  the  angle  of  scattering  of  fast  electron  in  inelastic  collisions 


Kig.9.  Average  values  of  cosine  of  the  angle  of  scattering  -  cos#*  and  of  direction  cosine  -  j.ik 

of  fast  electrons  in  e-beam  with  To=10  keV  depending  on  the  number  of  inelastic  collisions  of 
electrons  with  molecules  of  air. 
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Fig.  10,  The  average  energy  of  electrons  in  e-beam  with  To=10  keV  depending  on  the  number  of 
inelastic  collisions  of  electrons  with  molecules  of  air. 


x,  cm 

Fig.l  1 .  Spatial  distribution  of  energy  deposited  in  air  by  e-beam  with  T«=2keV,  No=3-10lr’cin  3. 
Red  curve  is  calculated  under  equation  (14),  black  curve  is  calculated  in  assuming  that  electron 
moving  direction  is  not  deflected  in  inelastic  collisions. 
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5.  Monte  Carlo  code  to  calculation  of  e-beam  propagation  in  nonuniform 
supersonic  flow  in  presence  of  magnetic  and  electric  fields  will  be  developed. 


In  the  stage  of  the  project  we  have  developed  the  Monte  Carlo  code  for  calculation  of  e- 
beam  propagation  in  nonuniform  supersonic  flow  in  presence  of  magnetic  and  electric  fields. 

Geometry  which  is  used  in  modeling  of  electrons  trajectories  is  shown  in  Fig.l.  The 
uniform  magnetic  field  B  directed  along  the  OZ  axis  and  the  uniform  electric  field  E  directed 
along  the  OY  axis  are  considered  in  modeling.  The  electron  beam  enters  the  gas  flow  at  the 
surface  z=0.  Moving  direction  of  the  e-beam  coincides  with  the  magnetic  field  direction.  The 
velocity  of  electrons  in  the  e-beam  significantly  exceeds  the  velocity  of  the  gas  flow,  thus  the 
flow  velocity  does  not  influence  on  the  e-beam  propagation  and  it  is  not  considering  in  the 
Monte  Carlo  modeling.  The  flow  and  the  e-beam,  uniform  in  OX  and  OY  directions,  are 
considered  in  the  modeling.  The  flow  parameters  in  OZ  direction  are  also  considered  as  uniform 
in  the  modeling  of  random  trajectories  of  electrons.  In  the  case  when  parameters  of  the  flow 
depend  on  the  z  coordinate  we  will  transform  results  of  the  Monte  Carlo  computations  by  using 
exact  relations  obtained  in  the  previous  stage.  The  Monte  Carlo  code  allows  us  to  calculate  the 
spatial  distribution  of  the  energy  deposited  by  e-beam  in  the  flow  and  any  characteristics  of 
electrons  in  the  e-beam. 

Trajectory  of  each  electron  in  the  Monte  Carlo  modeling  is  tracing  from  a  moment  of  its 
input  into  the  flow  with  initial  energy  To  up  to  a  moment  when  its  energy  becomes  less  than  the 
threshold  energy  Tsh.  Another  reason  to  stop  tracing  an  electron  is  implemented  when  the 
electron  gets  out  from  the  flow.  The  number  of  electrons  traced  in  the  Monte  Carlo  modeling  is 
called  as  the  number  of  histories.  Both  elastic  and  nonelastic  collisions  are  considered  in  the 
modeling. 

The  nonelastic  scattering  of  electrons  is  described  in  the  model  approach  which  was 
developed  in  the  third  stage  of  the  project.  The  differential  cross  section  describing  nonelastic 
losses  of  energy  of  fast  electrons  is  determined  in  this  case  by  the  equation: 


do.. 


Inr^mo 


0) 


dW  J32 

where  T  is  the  energy  of  electron,  W  is  the  energy  lost  by  electron  in  nonelastic  collision 
with  molecule  of  air,  ro  is  the  classical  electron  radius,  me  is  electron  mass,  c  is  velocity  of  light, 

t  =  t/ mec2  ,/3  =  v/c .  For  the  function  <f>(T,W )  the  fourth  approximation  is 


used:  (/>  =  <f>4(T,  IV)  = - — .  Detailed  description  of  the  functions  <f>4(T,  W )  and  r/4(T) 

(W  +  tj4(T))2 

is  (W  +  r]4(T))  given  in  previous  reports.  According  to  [1]  the  cosine  of  the  angle  of  scattering  of 
fast  electrons  in  nonelastic  collisions  with  molecules  can  be  calculated  under  the  formula: 

cos 0  ~  J(l-u')(r  +  2)/[(l-ir)r  +  2]  ,  (2) 

where  w=W/T  is  relative  value  of  the  energy  lost  by  electron  in  collision  with  molecule  of 
gas  (W  =  T-T),  T'  is  the  energy  of  electron  after  collision. 

The  differential  cross  section  of  elastic  scattering  of  fast  electrons  by  atoms  is  determined 
in  using  numerical  calculations  from  paper  [2],  We  have  proposed  new  expression  (3)  for  the 
differential  cross  section  which  provides  the  same  angle  dependency  as  one  traditionally  used  in 
approximate  calculations,  see  e.g.  [1,3].  The  newly  developed  expression  is  simpler  to  use  in 
approximating  results  of  numerical  calculations. 

doe,  _  °ei('1)  h(T){l+b( T)/2)  (3) 

dQ  2  n  (l-cos(#)  +  A(7’)): 

It  is  evident  that  equation  (3)  satisfies  the  known  dependency  of  the  total  cross  section  of 
elastic  scattering  of  electrons  oe/(T )  depending  on  the  electron  energy  T: 


43 


\^sLdCl  =  f^^-2;rsin(0V# 
J  dO.  {  dO.  w 


M7') 


The  function  b  in  equation  (3)  is  fitted  to  provide  coincidence  of  dependency  of  the  first 
transport  cross  section  cr^T)  calculated  in  [2]  with  results  of  calculation  the  function  <Ji(T)  in 
using  equation  (3)  for  the  differential  cross  section. 


<r,(p  =  J(l-<»(0))^<ra  =  ffo(7>»(7-) 
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(4) 


Thus  function  b(T)  can  be  determined  by  solving  equation  (4)  at  known  functions  of 
<Ji{T)  and  cjq(T)  .  In  order  to  calculate  the  total  cross  sections  and  the  first  transport  cross 


sections  for  elastic  scattering  of  fast  electrons  by  molecules  O2  and  N2  we  use  corresponding 
dependencies  for  atoms  O  and  N  which  are  presented  in  [2]  and  multiply  them  by  2.  This 
approach  is  frequently  used  in  Monte  Carlo  simulations  [1].  The  obtained  dependencies  of  total 
cross  sections  for  elastic  scattering  of  electrons  by  molecules  practically  coincide  with 
dependencies  used  in  paper  [3].  The  energy  losses  of  electrons  in  elastic  scattering  by  molecules 
are  significantly  less  than  ones  in  nonelastic  scattering  and  hence  they  are  not  included  in 
calculations. 

Random  magnitudes  of  energy  losses  W  of  electrons  in  an  individual  nonelastic  collision 
and  random  magnitudes  of  the  angle  of  scattering  of  electrons  in  an  individual  elastic  collision 
can  be  expressed  in  terms  of  random  number  $  uniformly  distributed  in  the  interval  [0,  1],  In  the 
Monte  Carlo  simulation  the  random  number  £  is  generated  by  a  random-number  generator. 
According  to  [4],  the  random  number  $  and  the  random  magnitude  of  any  variable  3  which  is 
changing  in  the  range  [«9min ,  <9max  ]  anc*  has  a  probability  density  p{3)  are  related  by  the 


equation:  J  p(3')d3'=  E, .  Thus  in  case  of  nonelastic  collision  the  function  can  be  determined  by 


solving  the  equation  (5): 


W(() 


J 

]^,(T,W’)dW 


=4 


(5) 


In  case  of  elastic  collision  of  electron  with  molecule  the  dependency  of  the  scattering 
angle  on  the  random  number  E,  can  be  determined  by  solving  equation  (6),  see  e.g.  [5]: 


m 


nil  1 


=  4 


(6) 


The  free  path  of  electrons  X  between  successive  collisions  with  air  molecules  has  a 
random  magnitude  too.  According  to  [4]  the  free  path  of  electrons  X  and  the  random  number  E, 
are  connected  by  the  equation  (7): 


A  = 


1 


■In 


(- 


(7) 


Average  free  path  of  electrons  in  this  case  is  equal  to  Ag  =  l/N(aej  +crn ) 

If  fast  electrons  propagate  in  the  gas  in  the  absence  of  electric  and  magnetic  fields  then 
trajectories  of  electrons  between  collisions  are  straight  lines.  In  conditions  of  MHD  controlled 
inlet  the  crossed  magnetic  and  electric  fields  are  imposed  on  the  flow.  Thus  in  general  case  the 
trajectory  of  electron  between  collisions  with  molecules  of  gas  differs  from  the  straight  line.  In 
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order  to  define  this  trajectory  it  is  necessary  to  solve  the  motion  equations  for  electrons  in 
electric  and  magnetic  fields.  It  is  necessary  to  note  that  high  energy  electrons  have  velocities 
which  are  comparable  with  velocity  of  light.  For  example,  electron  with  an  energy  T=100  keV  is 
characterized  by  the  ratio  of  its  velocity  to  the  velocity  of  light  p=v/c=0.548.  So,  in  general  case, 
the  relativistic  approach  needs  to  be  used  in  considering  an  electron  motion  between  collisions. 
According  to  [6]  the  vector  form  of  differential  equation  for  velocity  of  charged  particles  in 
electric  and  magnetic  fields  in  relativistic  approach  is  determined  by  equation  (8a): 


d\ 

~dt 


e 

m 


E  + vxB — --(v-E) 
c 


(8a) 


The  initial  condition  is  evident:  v(t  =  to)  =  Vo-  The  trajectory  of  electrons  is  determined  by 
the  motion  equation: 


dr 

dt 


v 


(8b) 


with  initial  condition  r(t  =  to)  =  ro. 

Fig.2  shows  trajectories  of  electrons  having  different  energies  at  the  initial  time  t  =  to  but 


the  initial  velocities  in  all  the  cases  have  the  same  direction: 


-M 


=  0.1, 


Vy(t  0) 


=  0.2 ,  initial 


v(to)  v(f0) 

coordinate  of  electrons  is  (0,0,0).  The  green  color  trajectory  in  the  Figure  corresponds  to 
conditions  with  B  =  0  and  E  =  0  .  One  can  see  that  trajectory  of  electron  significantly  depends  on 
the  energy  of  electron  and  on  the  magnetic  flux  density  magnitude.  Thus  in  conditions  of  MHD 
controlled  inlet  in  general  case  it  seems  that  influence  of  electric  and  magnetic  fields  on  the 
trajectory  of  electrons  between  collisions  need  to  be  taken  into  account.  Influence  of  this  effect 
on  the  spatial  distribution  of  energy  deposited  by  an  e-beam  in  the  flow  will  be  analyzed  in  the 
following  stage  of  the  project. 

Now  we  will  describe  the  modeling  algorithm  realized  in  the  Monte  Carlo  code 
developed  for  calculation  of  the  spatial  distribution  of  the  energy  deposited  by  e-beam  in  the  air 
flow  in  conditions  of  MHD  controlled  inlet.  At  the  initial  moment  at  the  point  where  an  electron 
enters  the  air  flow  the  initial  magnitudes  of  parameters  are  assigned  to  the  electron,  such  as  the 
initial  energy  T=To,  the  initial  coordinates  xo  ,yo,  zo  and  initial  direction  of  moving.  In  our 
investigations  we  assume  that  parameters  both  of  the  flow  and  the  e-beam  are  uniform  in  OX  and 
0Y  directions.  So  the  spatial  distribution  of  the  energy  deposited  by  the  e-beam  will  be  a  function 
only  z  coordinate.  It  is  evident  that  in  our  consideration  initial  values  of  coordinates  x  and  y  are 
of  no  importance  in  calculations  and  we  will  assume  xo=0  ,  yo=0.  The  surface  from  which 
electrons  enter  the  flow  corresponds  to  z=0,  so  the  initial  value  for  z  coordinate  of  electron  is 
zo=0.  The  value  of  cosine  of  the  angle  between  the  initial  moving  direction  of  an  electron  and  the 
0Z  axis  is  |J  o=l. 

The  trajectory  of  an  electron  between  two  successive  collisions  is  defining  from  solution 
of  the  set  of  equations  (8a,  8b).  When  J3  «  1  the  solution  can  be  determined  in  analytical  form, 
otherwise  the  equations  need  to  be  solved  numerically.  The  time  of  an  electron  motion  between 
the  successive  collisions  can  be  approximately  determined  by  the  relation :t  =  A/v,  where  v  is 
velocity  of  electron  after  the  preceding  collision  with  a  molecule  (estimations  shows  that  in 
conditions  of  MHD  controlled  inlet  the  modification  of  absolute  value  of  the  electron  velocity 
along  its  free  path  is  negligible).  In  order  to  determine  a  random  magnitude  of  the  free  path  of  an 
electron  X  the  random  number  £  is  generating  and  then  the  equation  (7)  is  used.  The  calculated 
trajectory  of  the  electron  allows  one  to  define  coordinates  of  the  point  in  which  the  electron  will 
be  scattered  by  a  molecule.  The  moving  direction  of  the  electron  at  the  moment  of  collision  will 
be  determined  too.  After  that  the  next  random  number  £  is  generated  in  order  to  define  what 


type  of  collision  occurs.  If  the  inequality  £  <  cre/  /(cre/  +<7n)  is  true  the  collision  is  elastic 
otherwise  the  collision  is  nonelastic.  A  next  random  number  t;  is  generated  to  model  a  collision. 
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In  case  of  elastic  collision  the  scattering  angle  6  is  determined  by  substituting  the 
random  number  t;  in  equation  (6).  The  energy  W  lost  by  the  electron  in  elastic  collision  is 
considered  to  be  zero,  W=0.  In  case  of  nonelastic  collision  the  magnitude  of  the  energy  lost  by 
the  electron  W  is  determined  under  equation  (5).  The  scattering  angle  0  of  electron  in  the 
nonelastic  collision  is  determined  by  equation  (2).  The  azimuth  angle  of  scattering  y/  is 
calculated  both  for  elastic  and  nonelastic  collision  by  using  a  next  generated  random  number  £ 
under  the  relation:  y/  =  2n<; .  The  energy  of  electron  after  a  collision  T’  is  determined  by  the 
equation  V  -  T-W.  After  the  individual  collision  of  the  electron  has  been  modeled  the  moving 
direction  for  the  electron  changed  in  the  collision  will  be  determined  by  using  standard 
geometrical  transformations  [1],  These  modified  parameters  of  the  electron  will  be  considered  as 
initial  parameters  in  modeling  the  electron  trajectory  up  to  subsequent  collision  with  a  molecule. 
This  process  of  modeling  is  repeated  many  times.  The  modeling  of  trajectories  of  an  individual 
electron  will  be  over  in  two  cases:  the  electron  energy  becomes  less  than  the  threshold  energy 
T’  <  T,i„  or  the  electron  leave  the  investigated  region  under  condition  z<0. 

In  order  to  verify  the  developed  Monte  Carlo  code,  we  compared  the  results  of  Monte 
Carlo  computations  with  exact  results  obtained  analytically  in  previous  stages  of  the  project.  The 
features  which  were  not  taken  into  account  in  obtaining  the  analytical  results  were  not  taken  into 
account  in  the  Monte  Carlo  computations  too.  Fig.3  demonstrates  dependency  of  average  energy 
of  electron  depending  on  the  number  nonelastic  collisions  with  molecules  of  air.  Figs.4-6 
demonstrate  the  energy  distribution  function  of  electrons  after  k  nonelastic  collisions  with 
molecules  of  air.  Fig.7  demonstrates  the  spatial  distribution  of  the  energy  deposited  by  the  e- 
beam,  calculated  in  approximation  u=l.  Only  nonelastic  collisions  were  considered  in  this 
approach.  Fig.8  demonstrates  the  average  value  of  the  cosine  of  the  scattering  angle  of  an 
electron  in  nonelastic  collisions  depending  on  the  number  of  collisions  accomplished  by  the 
electron.  One  can  see  that  Monte  Carlo  computations  are  in  a  good  accordance  with  the  exact 
analytical  results.  Thus  we  can  conclude  that  the  developed  Monte  Carlo  code  is  useable  and  it 
can  be  used  for  verifying  approximate  analytical  solution  obtained  in  the  previous  stage.  Fig.9 
demonstrates  results  of  calculations  obtained  in  consideration  only  nonelastic  collisions  of 
electrons  with  molecules  of  air.  In  the  approach  1  the  direct-forward  motion  of  e-beam  is 
considered  (ju  =  1).  In  the  approach  2  the  angle  scattering  of  electrons  in  nonelastic  collisions  are 
considered.  One  can  see  that  analytical  solution  in  the  approach  2  significantly  differs  from  the 
Monte  Carlo  computation.  So  we  can  conclude  that  some  of  approximations  which  were  used  in 
obtaining  the  analytical  solution  for  the  approach  2  need  to  be  corrected.  The  reasons  which  are 
responsible  for  such  difference  will  be  analyzed  in  the  following  stage  of  the  project.  Figs.  10-12 
demonstrate  comparison  of  experimentally  measured  in  [7,  8]  the  spatial  distributions  of  energy 
deposited  by  electrons  in  the  air  and  in  the  nitrogen  with  our  Monte  Carlo  computations.  One  can 
see  very  good  accordance  between  numerical  and  experimental  results.  Thus  we  can  conclude 
that  the  developed  Monte  Carlo  code  is  adequately  describes  the  process  of  energy  deposition  by 
e-beam  and  it  can  be  used  to  calculate  the  energy  deposition  profiles  in  conditions  of  MHD 
controlled  inlet. 
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CONCLUSIONS 

In  the  fifth  stage  of  the  project  the  Monte  Carlo  code  was  developed  for  calculation  of  e- 
beam  propagation  in  the  air  flow  in  conditions  typical  for  the  MHD  controlled  inlet.  This  code 
was  verified  by  comparison  of  results  of  the  Monte  Carlo  computations  both  with  exact 
analytical  results  obtained  in  previous  stages  of  the  project  and  with  known  experimental  data  of 
spatial  distributions  of  the  energy  deposited  by  electrons  in  the  air  and  nitrogen.  Good 
accordance  between  the  compared  results  demonstrates  that  the  developed  Monte  Carlo  code  is 
adequately  describes  the  process  of  e-beam  propagation  in  the  flow  and  it  can  be  used  to 
calculate  the  spatial  distribution  of  the  energy  deposited  by  e-beam  in  the  flow  in  conditions  of 
MHD  controlled  inlet. 
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Electron  beam 


Fig.].  The  geometry  considered  in  modeling  the  e-beam  propagation  in  conditions  of  MHD 
controlled  inlet 


p=0.  i 
|3=0. 5 


0.3  -0. 1 

Fig.2  Trajectories  of  high-velocity  electrons  in  magnetic  and  electric  fields  at  various  values  of 
P=v/c.  B=BZ=3T,  E=Ey  =  30  V/cm,  vr(/0)/v(/0)-0  1,  vF(/0)/v(/0)  =  0.2,  N=10!W3. 

The  green  curves  correspond  to  B=0,  E=0. 
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6.  Comparative  analysis  of  space  distributions  of  power  density  spent  on  flow 
ionization  by  e-beam  calculated  both  in  Monte  Carlo  code  and  in  Boltzmann 
kinetic  equation  solution.  Improvement  of  mathematical  model  of  e-beam 
propagation  in  nonuniform  supersonic  flow  will  be  done.  Rapid  method  to 
calculate  space  distribution  of  electric  conductivity  at  given  parameters  of  e-beam 
and  structure  of  gas-dynamic  flow  in  inlet  will  be  developed. 

In  the  stage  of  the  project  the  comparative  analysis  of  spatial  distribution  profiles  of  the 
energy  deposited  by  e-beam  in  the  air  flow  calculated  both  in  Monte  Carlo  code  and  in  analytic 
solution  of  the  Boltzmann  kinetic  equation  was  made.  Improvement  of  analytical  modeling  was 
made  by  developing  the  method  of  transformation  of  analytical  dependencies.  Rapid  method  for 
calculation  of  spatial  distribution  of  power  density  deposited  by  e-beam  in  the  air  flow  and  hence 
for  calculation  of  electron  concentration  and  electric  conductivity  in  a  nonuniform  supersonic 
flow  at  given  parameters  of  e-beam  is  developed. 

The  spatial  distribution  profiles  of  the  energy  deposited  by  e-beam  in  the  air  flow  have 
been  calculated  in  using  Monte  Carlo  code  developed  in  the  previous  stages  of  the  project. 
Calculations  were  made  for  set  of  initial  energies  of  electrons  in  the  e-beam  in  the  range  from  2 
to  100  keV.  Comparison  of  results  of  the  Monte  Carlo  calculations  with  results  of  approximate 
analytical  solution  in  which  the  deflection  of  electrons  in  scattering  is  taken  into  account  has 
been  made.  It  is  shown  that  the  main  source  of  discrepancy  between  the  Monte  Carlo 
calculations  and  the  approximate  analytical  solution  is  the  approximation 

g(T',ju'—. >  T,ju)  =  g(T'^>  T)-g(ju'—. >  //)  which  was  used  to  account  analytically  the  angle 
scattering  of  electrons  in  inelastic  collisions  with  molecules  of  air.  This  inference  is  illustrated  in 
Fig.l  where  dependencies  of  the  average  value  ///(  of  the  direction  cosine  of  the  angle  of  fast 
electrons  in  the  e-beam  after  k  inelastic  collisions  are  shown  as  a  function  of  k.  One  can  see  that 
calculated  in  the  Monte  Carlo  code  (the  curve  1)  differs  sharply  from  approximate  analytical 

results  (the  curve  3).  The  curve  2  shown  in  the  figure  is  calculated  by  Monte  Carlo  code  in  using 
approximations  which  adequate  to  the  approximations  used  in  obtaining  analytical  results. 
Namely,  in  order  to  calculate  the  dependency  2  in  the  Monte  Carlo  code  the  approximation  was 
included  that  the  energy  loss  and  the  scattering  angle  in  electron-molecule  collisions  are 
independent.  One  can  see  that  analytical  results  and  Monte  Carlo  calculations  obtained  in  similar 
approaches  are  practically  coincides.  But  the  approximation 
g(T' T,ju)  #jg(r'— »  T)-  g(ju'—>  ju)  which  was  used  in  analytical  calculations  is  not 
adequate  for  describing  multiple  collisions  despite  the  fact  that  the  energy  loss  and  the  angle  of 
scattering  of  electron  in  individual  collisions  are  usually  very  small.  In  order  to  obtain  more 
reliable  method  of  analytical  calculation  of  profiles  of  the  energy  deposition  we  have  used  the 
Monte  Carlo  computations  to  correct  the  developed  analytical  approaches. 

Comparison  of  the  Monte  Carlo  computations  and  exact  analytical  solution  obtained  in 
the  “straightforward”  approach  (with  //=  1)  demonstrates  correlation  between  the  shapes  of 
dependencies  of  spatial  distribution  of  energy  deposited  by  e-beam  for  various  values  of  initial 
energy  of  electrons  in  the  e-beam.  Two  examples  of  such  comparisons  are  shown  in  Figs.2,  3  for 

Td= 4keV  in  Fig.2  and  7o=80  keV  in  Fig.3.  Here  z  is  dimensionless 

X 

coordinate:  z  =  E  •  jTVfx'jJx',  where  N(x)  is  concentration  of  molecules  depending  on  the  spatial 
0 

coordinate  x,  E  is  the  total  cross-section  of  inelastic  scattering  of  fast  electrons  on  molecules  of 
air  calculated  in  approach  of  the  modeling  differential  cross-section  developed  in  the  third  stage 

of  the  project.  In  case  of  constant  concentration  of  molecules  z  =  NXx.  In  order  to  make  a 
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correction  of  analytical  results  it  is  more  convenient  to  compare  analytical  and  Monte  Carlo 
results  not  for  dT  /  dz  dependencies  but  for  the  energy  flux  T(z)  which  is  calculated  in  terms  of 

T  1 

the  flux  density  d>(z,//,T )  by  the  relation:  T(z)  =  ^TdT  j  juQ>(z,jU,T)djU .  In  the  analytical 

0  -1 

solution  for  the  straightforward  approach,  obtained  in  the  third  stage  of  the  project,  the  energy 

T 
J  0 

flux  can  be  calculated  under  equation T(z)  =  ^X/C(z)  jV •  F;c(T)dT  =^X !c(z)Tj. ,  where  7)c  is 

k  o  k 

average  energy  of  electron  after  k  inelastic  k  collisions,  X^{z)  =  z1'  exp(-z)/ k\ .  In  order  to 

obtain  the  energy  flux  T(z)  by  Monte  Carlo  method  we  are  integrating  the  dependencies  of  dT/dz 

00 

which  are  obtained  in  the  Monte  Carlo  calculations:  T(z)  =  ^ (dT / dz')dz' .  The  energy 

0 


1 

fluxes  transferred  by  an  e-beam  (the  initial  flux  is  assumed  to  be  unit  i.e.  J  juO(0,jU,T)djU  =  1 ) 

-1 

are  shown  in  Figs.4,5  for  two  values  of  initial  energies  depending  on  the  distance  z  the  e-beam  is 
passing  through.  Red  curves  correspond  to  analytical  results  and  black  curves  correspond  to 
Monte  Carlo  calculations.  One  can  see  that  the  deflection  of  electrons  in  elastic  and  inelastic 
scattering  taken  into  account  in  the  Monte  Carlo  calculations  leads  to  faster  decreasing  of  the 
energy  flux  in  comparing  with  the  straightforward  analytical  approach  where  motion  direction  of 
electrons  does  not  change  in  collisions.  The  energy  flux  obtained  in  the  Monte  Carlo  calculations 
at  z=0  is  also  less  than  one  in  the  analytical  approach.  Such  difference  results  from  the  fact  that 
part  of  electrons  from  e-beam  is  returned  to  the  surface  z=0  when  their  energies  are  high  enough. 
Thus  a  part  of  the  energy  flux  of  e-beam  leaves  the  air  flow  together  with  the  returned  electrons. 

In  order  to  improve  analytical  modeling  we  have  developed  method  of  transformation  of 
analytical  dependencies  obtained  in  the  straightforward  approach.  The  dependencies  of  energy 
fluxes  obtained  analytically  in  the  straightforward  approach  -  TSF(z)  were  compared  with 
corresponding  dependencies  obtained  in  the  Monte  Carlo  calculations  -  TMC(z)  for  several  values 
of  initial  energy  of  electrons  To  in  the  e-beam.  In  result  of  the  comparison  the  equations  were 
obtained  which  allows  us  to  calculate  analytically  dependencies  of  T(z)  and  dT/dz  which  are  in 
good  accordance  with  the  Monte  Carlo  calculations.  Basic  features  of  such  comparison  are 
shown  in  Fig.6.  We  make  the  coordinate  transformation,  by  solving  the  equation 
Tsf(z  (z))  =  Tmc(z).  The  parameter  z0,  shown  in  Fig.6.,  is  a  solution  of  the  equation 

Tsf(z0)  =  Tmc(0).  In  result  of  calculations  at  several  values  of  To  changing  in  the  range  of 


2M00keV  we  have  obtained  approximating  function  for  z0 : 

2- 105  •  Tv 2 

l  = _ -fq _  (i) 

0  2.659- 104  -  39.12  -7^2  +  0.2869 -T^2 

Parameter  To  in  Eq(l)  is  measured  in  keV.  The  equation  (1)  can  be  used  for  To  changing 
in  the  range  2<T0<  100.  The  function  zq(Tq)  is  shown  in  Fig.7;  dots  in  the  figure  correspond  to 

solution  of  Tsf  (  zq  )  —  TMC  (0)  at  given  energies.  According  to  Fig.8  the  dependencies  of 


z/zQon  the  z/zq  ratio  are  practically  coincides  for  various  initial  energies  of  electrons  in  the 

investigated  range  2^-100  keV.  The  approximating  function  (Eq(2)),  shown  in  the  Figure,  was 
obtained  to  describe  this  dependency. 
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Equations  (1,2)  allows  us  to  obtain  analytical  dependencies  of  the  spatial  distribution  of 
the  energy  deposited  by  e-beam  dT  /  dz  in  the  flow  by  transforming  the  dependencies  obtained  in 
the  straightforward  approach.  Corresponding  relation  for  this  transformation  is  presented  below: 

dJ,  ■.  £/Tsf  ^  dz  ... 

(3) 

Figs. 9, 10  demonstrate  that  the  spatial  distributions  calculated  analytically  by  using  the 
transformation  method  are  in  good  accordance  with  results  of  Monte  Carlo  calculations.  The 
next  step  was  made  in  the  stage  of  the  project  to  simplify  equations  for  calculation  of  average 
energy  of  electron  after  k  inelastic  collisions  Tk,  in  order  to  make  more  rapid  method  for 

calculation  of  dT/dz.  In  the  third  stage  of  the  project  the  average  energy  7)c  was  calculated  by 
using  the  following  relations: 

7.  -  j  ?*•/*. W/ 

r.  °  (4) 

F, (T) = Iff (r -> (ryr,  f,(t) =s(t-t„) 


F, (?) = Jff (r -> (T)dT\  F„(T) =S(T-T„) 

T 

So,  in  order  to  find  average  energy  7^  of  electron  after  k  collisions  we  need  to  find 

previously  the  energy  distribution  function  Fk(T)  of  electron  after  k  collisions.  In  this  stage  of  the 
project  we  have  shown  that  the  average  energy  can  be  calculated  without  calculation  of  the 
energy  distribution  function.  We  introduce  the  function  7).  (T)by  the  equations: 

w)=\ur)-x(T^ryr  at  k>] 

T  (5) 

7;(7’)=Jr-g(r^ryr 

It  is  shown  that  the  average  energy  of  electrons  having  initial  energy  To  can  be 
determined  in  terms  of  introduced  function  by  the  relation:  7)c  =  7)c  (Tq  ) .  Some  dependencies  of 

Tj((T)avQ  shown  in  Fig.l  1.  One  can  see  that  Tk(T)  is  monotone  increasing  function  in  contrast  to 
the  nonmonotonic  energy  distribution  function  Fk(T)  and  it  is  easy  to  be  used  for  calculations  of 
TSF  (z)  by  the  equation:  TSF  (z)  =YJXk(z)Tk(T0) . 

k 

The  following  relation  determines  algorithm  of  transformation  of  obtained  results  to  real 
coordinate,  instead  of  dimensionless  z  in  case  of  e-beam  propagating  along  the  y  coordinate: 


{y)  =  ~T-  z  =  Z-\N{y’)dy'  -  — r  =  -N(y)Z 


If  the  spatial  distribution  of  a  molecular  concentration  along  the  e-beam  path  N(y)  is 
known  we  can  calculate  the  spatial  distribution  of  the  power  density  deposited  by  e-beam  in  the 
flow  by  the  relation:  qion  =  (//,  /  e)dT  /  dy ,  where  jb  is  the  e-beam  current  density,  e  is  electron 
charge. 

In  considering  nonuniform  flow  the  flow  parameters  usually  will  be  depending  not  only 
on  the  y  coordinate.  In  2D  approach  we  will  have  dependency  N(x,y).  Possible  configuration  of 
the  region  of  the  e-beam  power  deposition  in  nonuniform  flow  with  oblique  shock  is  shown  in 
Fig.  12.  Here  we  consider  a  uniform  magnetic  field  directed  along  the  y  axis  and  the  e-beam 
propagating  along  the  magnetic  field,  the  width  of  the  e-beam  is  xj-  -  vy, .  In  this  case  the  power 

density  deposited  by  e-beam  in  the  flow  will  be  a  function  of  x  and  y  coordinates.  And  it  can  be 
calculated  by  the  relation: 
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Results  of  calculations  corresponding  to  supersonic  flow  around  the  wedge  with  angle  of 
6.5°  are  shown  in  Figs.  13-15.  Calculations  were  made  at  free  stream  dynamic  pressure  gx=751cPa 
with  Mach  number  of  incoming  flow  M^=6.  Figs.  13- 15  demonstrate  correspondingly  the  spatial 
distribution  of  the  power  density  qion  deposited  by  e-beam  in  the  flow,  the  electron  concentration 
ne  sustained  by  the  e-beam,  and  the  flow  conductivity  CT.  The  influence  of  effective  electric  field 
on  concentration  of  electrons  and  the  flow  conductivity  was  taken  into  account  in  accordance 
with  results  presented  in  previous  stages  of  the  project. 


CONCLUSIONS 


In  this  stage  of  the  project  the  comparative  analysis  of  spatial  distribution  profiles  of  the 
energy  deposited  by  e-beam  in  the  air  flow  calculated  both  in  Monte  Carlo  code  and  in  analytic 
solution  of  the  Boltzmann  kinetic  equation  was  made.  The  method  of  transformation  of 
analytical  dependencies  obtained  in  the  straightforward  approach  is  developed.  Rapid  method  for 
calculation  of  spatial  distribution  of  power  density  deposited  by  e-beam  in  the  air  flow  and  hence 
for  calculation  of  electron  concentration  and  electric  conductivity  in  a  nonuniform  supersonic 
flow  at  given  parameters  of  e-beam  is  developed.  Calculations  of  power  density  deposited  by  e- 
beam  in  the  air  flow,  the  electron  concentration  and  conductivity  of  the  flow  were  made  in 
conditions  of  nonuniform  flow  with  oblique  shock. 


Fig.  1 .  Average  value  of  direction  cosine  of  fast  electrons  after  k  inelastic  collisions  at  Trj=  2  keV. 

I  -  Exact  Monte  Carlo  calculations;  2  -  Approximate  analytical  solution;  3  -  Approximate 
Monte  Carlo  calculations 
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Fig.2  Spatial  distribution  of  the  energy  deposited  by  e-beam  with  7^=4  keV.  Black  curve  is 
Monte  Carlo  calculations;  red  curve  is  analytical  calculations  in  the  “straightforward”  approach 


Fig.3  Spatial  distribution  of  the  energy  deposited  by  e-beam  with  7^=80  keV.  Black  curve  is 
Monte  Carlo  calculations;  red  curve  is  analytical  calculations  in  the  “straightforward”  approach 
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Fig.4.  Energy  flux  transferring  by  e-beam  at  Trj=4  keV.  Black  curve  is  Monte  Carlo  calculations; 
red  curve  is  analytical  calculations  in  the  “straightforward”  approach 


Fig. 5.  Energy  flux  transferring  by  e-beam  at  7^=100  keV.  Black  curve  is  Monte  Carlo 
calculations;  red  curve  is  analytical  calculations  in  the  “straightforward”  approach 
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Fig.6.  Algorithm  of  transformation  of  analytical  dependency 


Fig.7.  Approximating  function  z0(7^) 
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Fig. 8.  Approximating  function  for  z/z0  ratios  obtained  for  various  values  of  To. 


Fig.  9.  The  spatial  distribution  of  the  energy  deposited  by  e-beam  calculated  analytically  by  using 
the  transformation  method  (black  curve)  at  7^10  keV.  Red  curve  is  Monte  Carlo  calculation. 
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Fig.  10,  The  spatial  distribution  of  the  energy  deposited  by  e-beam  calculated  analytically  by 
using  the  transformation  method  (black  curve)  at  T(f=40  keV.  Red  curve  is  Monte  Carlo 
calculation. 


Fig.  1  la.  The  function  Tk  (7')  at  various  numbers  of  collisions  k. 
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Fig.  1  lb.  The  function  7’  (/)  at  various  numbers  of  collisions  k. 


Fig.  12.  Region  of  the  e-beam  power  deposition  in  nonuniform  flow  with  oblique  shock. 
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Fig  13  Power  density  deposited  by  e-beam  in 
the  air  flow  at/&=10  mA/cm2,  To=40  keV, 
B=2J 


Fig.  14.  Electron  concentration  in  the  air  flow 
sustained  by  e-beatn  at  //,= 1 0  mA/cm2,  7o=40 
keV,  B= 2T 


Fig.  15.  Conductivity  in  the  air  flow  sustained  by  e-beam  at y&=l  0  m A/cm 2,  7>~40  keV,  B= 2T 
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7.  Mathematical  model  of  MHD-controlled  inlet  with  nonequilibrium  MHD 
generator  using  e-beam,  as  ionizer,  will  be  developed  in  2D  Euler  approach. 
Iterative  scheme  will  be  used,  which  takes  into  account  that  both  flow  structure 
influence  on  conductivity  distribution  and  nonequilibrium  MHD  generator  modify 
flow  structure. 

“Cold  air  flow”  as  used  herein  for  the  purposes  of  this  investigation  implies  the  flow  with 
thermal  ionization  degree  insufficient  for  realization  of  MHD  control  principles.  The  importance 
of  this  investigation  stems  from  the  fact  that  in  case  of  insignificant  natural  flow  conductivity  it 
is  necessary  to  take  measures  to  increase  its  ionization  degree.  The  estimates  show  that  such  state 
of  affairs  is  typical  for  hypersonic  flight  conditions,  i.e.  ionizer  is  a  crucial  part  of  any  MHD 
system. 

The  energy  required  for  production  of  one  electron-ion  pair  we  shall  call  “cost  of 
ionization”  and  designate  W,.  Analysis  of  the  known  ionization  techniques  show  that  from  the 
energy  point  of  view  an  electron  beam  is  the  least  expensive  means  for  gas  ionization.  A  price  of 
ionization  in  this  method,  with  electron  energy  in  a  beam  exceeding  1  KeV  makes  W)  ~  34  eV. 
This  value  exceeds  insignificantly  the  ionization  energy  for  air  molecules.  High-voltage  impulse 
discharge  has  larger  ionization  price  making  around  66  eV  for  air. 

Irrespective  of  ionization  technique,  we  shall  characterize  ionizer  operation  by  power  q, 
deposited  into  unitary  flow  volume.  Then  ionization  velocity  in  the  unitary  volume  will  be 

defined  by  j.  =  —  (1) 

Wi 

For  example,  at  the  gas  ionization  by  the  electron  beam,  value  q,  is  defined  by  the  ratio: 

qi-^D(Eb)P , 

vj 

where  j  b  is  current  density  in  the  beam,  D(E  b)  -  electron  breaking  ability,  E  b  -  electron  energy 
in  a  beam. 

Ionized  gas  conductance  is  determined  by  electron  concentration  that  at  stationary  state  and 
preset  velocity  is  gauged  through  the  velocity  of  their  destruction  processes.  In  the  case  that  the 
two-part  dissociative  recombination  (e.g.  with  positively  charged  molecular  ions)  is  a  prevailing 
process  of  electron  destruction,  steady-state  value  of  electron  concentration  is  specified  by 


ne 


g/ 

\  kdr  Wj  ’ 


(2) 


where  kdr  is  two-part  electron-ionic  recombination  rate  constant. 

As  a  case  in  point,  Fig.l  provides  dependence  of  electron  concentration  on  specific  costs 
for  ionization  qi  calculated  in  the  event  of  electron  beam  usage  as  a  means  of  ionization. 
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q W/cm3 

Fig.l.  Dependence  of  electron  concentration  in 
air  plasma  on  energy  deposited  through 
electron  beam  into  unitary  volume  in  a  unit  of 
time 


Fig.2.  Dependence  of  air  plasma  ionization 
degree  on  energy  deposited  through  electron 
beam  into  unitary  volume  in  a  unit  of  time 


The  presented  dependence  is  well  approximated  by  the  expression 

ne  =  1.124  l012-^Q/>  cm'3, 


(3) 


where  q,  is  given  in  W/cm  .  Approximation  error  here  does  not  exceed  30  %. 

It  follows  from  (1)  that  flow  ionization  degree  can  be  presented  in  the  following  form: 


ne_  QjJjh 
no  V  kdr  Wj 


(4) 


2 

Computational  results  of  flow  ionization  degree  depending  on  the  value  q,/no  are  given  in  Fig.2. 
Numerical  computation  results  are  well  approximated  by  the  expression  below 


He 

n0 


=  1.17  io-5(l03Vn2oy 


2x0.6 


(5) 


'J  'J 

where  q,  is  given  in_W/cm  units  and  concentration  n0  -  in  cm'  .  Approximation  error  in  this  case 
does  not  exceed  9%. 

Formulae  (3)  and  (5)  can  be  used  for  prompt  evaluation  of  plasma  ionization  degree  and  electron 
concentration,  as  well  as  for  analytical  computation  in  the  sphere  of  parameter  values  meeting 
the  inequality: 

2-  10  2  <  (q;  /  n02)  •  1034  <20  (6) 

Plasma  conductivity  is  connected  with  electron  concentration  by  the  ratio 

o  =  e2/ie/meno/cc  ,  (7) 


where  e  is  electron  discharge;  me  -  electron  mass;  kc  -  electron  dissipation  rate  constant. 
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The  expression  is  deduced  from  (2)  and  (7)  that  bounds  flow  conductivity  to  specific  power 

— (8) 

menokc  v  k*  Wi 


spent  for  ionization:  CT  = 


It  is  evident,  that  MHD  generator  with  non-equilibrium  conductivity  could  be  realized  only 
in  the  event  when  the  energy  spent  for  flow  ionization  to  produce  required  conductivity  level 
would  not  exceed  energy  produced  by  MHD  generator.  Such  operation  conditions  for  the  MHD 
generator  with  ionizer  we  call  self-sustained  mode. 

Power  produced  by  MHD  generator  depends  on  flow  conductivity,  flow  gasdynamic  parameters, 
MHD  performance  and  MHD  generator  type.  Thus,  conditions  wherein  the  self-sustained  mode 
could  be  realized  depend  on  many  factors,  such  as:  MHD  generator  arrangement  within 
hypersonic  flight  vehicle  (HFV),  ramjet  geometry,  flight  trajectory,  MHD  performance,  ionizer 
type,  etc.  In  this  paper  the  version  is  looked  at  where  MHD  generator  of  the  ramjet  is  located 
upstream  of  the  combustor  and  MHD  accelerator  is  placed  downstream  of  the  combustion 
chamber. 

Let  us  identify  variation  region  of  MPCE  parameters  wherein  self-sustained  mode  for 
ionizer  and  MHD  generator  exists.  To  do  that,  we  need  to  correlate  power  spent  for  ionization 
with  energy  produced  by  MHD  generator.  Volume  density  of  the  power  produced  by  MHD 
generator  is  given  by  expression 

qg=k(l-k)aBV  ,  (9) 


where  k  is  load  factor  of  MHD  generator;  a  -  current  conductivity;  v  -  flow  velocity;  B  - 
magnetic  field  induction. 

As  flow  conductivity  depends  on  power  spent  for  flow  ionization,  specific  energy  produced 
by  MHD  generator  is  function  of  parameter  qi: 

=  -Sl-  (10) 

meriokc  V  kdr  Wi 


Some  results  of  numerical  computation  of  dependencies  of  MPCE  power  characteristics  on 
ionizer  parameters  and  its  location  points  for  various  parameters  of  the  free  air  flow  are 
presented  in  Figs.  3  -  5. 

Dependence  of  qg  -q,  (parameter  difference  of  power  produced  and  energy  deposition  into 
ionization)  on  parameters  q;  and  B  is  provided  in  Fig.3.  It  is  obvious,  that  there  is  a  limiting 
value  for  magnetic  induction  B  =  Bcr  for  which  power  produced  by  MHD  generator  exceeds 
ionization  cost:  qg  >  q,.  With  energy  deposition  increase,  critical  value  of  magnetic  induction 
increases.  Magnetic  induction  enhancement  at  the  fixed  ionization  cost  leads  to  the  increase  of 
power  produced  by  MHD  generator  and  expands  the  range  of  self-sustained  mode  of  MHD  with 
non-equilibrium  conductivity. 
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Fig.3.  Power  density  produced  by  MHD 
generator  less  power  density  spent  by  ionizer  (qg 
-  qi )  depending  on  parameters  qi_and  B.  MHD 
generator  is  placed  at  cross-section  2  ,  Mx.  =  6 


Fig.4.  Power  density  produced  by  MHD 
generator  with  non-equilibrium 
conductivity  less  flow  ionization  cost  at 
various  values  of  magnetic  induction  as 
function  of  power  deposited  into 
ionization. 


Dependence  of  qg  -  qi  is  not  monotonic.  There  exists  some  optimal  value  q,  wherein  value 
qg  -  q,  has  the  maximum  for  the  given  magnetic  induction  value  (Fig.4).  In  the  event  of  the 
expression  (8)  possible  usage  to  compute  ionized  flow  conductivity,  magnetic  induction  critical 
value  could  be  obtained  from  equation  (10): 


Bcr 


menokcVq'/kdrl/l// 

k(1-k)eV 


(11) 


ft  follows  from  the  presented  ratio,  that  critical  value  of  magnetic  induction  depends  on 
flow  parameters,  MHD  generator  and  ionizer,  as  well  as  on  velocity  constants  of  elementary 
plasma  processes.  On  the  other  hand,  at  magnetic  induction  fixed  value,  it  is  possible  to 
determine  power  density  critical  value,  wherein  specific  ionization  cost  are  equal  to  specific 
power  produced  by  MHD  generator  (qg  =  q,)-  From  expression  (10)  for  critical  value  of  specific 
power. 


qcr 


i 


kdr  Wi 


/c(i-/c)eV - 

V  men0k  c) 


(12) 


In  this  case  self-sustained  mode  for  MHD  generator  and  ionizer  is  realized  at  qcr  >  q,  >  0. 

From  the  relation  (12)  it  follows  that  the  existence  range  of  MPCE  self-sustained  mode 
expands  with  flow  velocity  increase  and  gas  concentration  decrease.  Thus,  the  faster  and  higher 
HFY  flies,  the  more  efficient  MHD  systems  application  on  its  board  is.  Magnetic  induction 
increase  also  leads  to  the  expansion  of  existence  range  for  self-sustained  mode  of  MPCE  and 
ionizer. 

Let  us  evaluate  the  value  qcr,  constant  values  for  elementary  processes:  kdr  =  1.5  10'  cm/s, 
kc  =2  10"8cm3/s,  W;=34  eV  for  the  flow  with  the  following  parameters:  v  =  2000  v/sec;  T  =  300K 
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at  values  B  =  2T,  k  =  0.5.  We  obtain  qcr  =  6.51  W/cm  .  Value  qcr  is  one  of  the  most  important 
characteristics  of  MPCE  operation  conditions. 

By  way  of  example,  shown  in  Fig.  5  is  effect  of  relation  qt  /qcr  on  density  distribution  at 
MHD-controlled  inlet  in  cross-section  of  the  channel  upstream  of  MPCE  combustor.  Values  q-, 
/qcr  are  indicated  adjacent  to  the  curves.  Flow  density  increases  with  q,  /qcr  rise.  Under  given 
conditions  for  this  computation,  significant  density  increase  is  achieved  even  at  q-,  /  qcr  =0.4. 
Further  increase  of  qt  /_qa-  relation  under  these  conditions  does  not  result  in  significant  variation 
of  flow  density  distribution.  Excessive  power  qg  -  qi  left  after  expenditures  for  ionizer  operation 
could  be  realized  on  board  HFV  for  other  needs. 


y.  arb.u.  W/CrTI 


Fig.  5.  Density  distribution  at  the  channel  cross- 
section  upstream  of  MPCE  combustor;  Mx  =6; 
dynamic  pressure  makes  3.7  •  104  Pa;  values  q, 
/qcr  are  indicated  adjacent  to  the  curves. 


Fig.6.  Region  of  parameters  B  and  q,  variation 
wherein  self-sustained  mode  for  MHD 
generator  with  non-equilibrium  conductivity 
and  required  parameter  value  of  MHD 
interaction  are  simultaneously  ensured.  MHD 
generator  is  located  at  the  sectionl.  Mx  =  6. 


Table  1  provides  values  qmaXG  =  qg-qi  which  is  maximum  power  taken  off  MHD  generator 
volume  unit  less  power  consumed  by  ionizer  at  different  values  of  temperature,  flow  velocity, 
magnetic  induction  magnitude  and  load  factor.  Working  gas  pressure  is  equal  to  0.1  atm. 
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Table  1 


Characteristics  of  the  specific  excessive  power  of  MHD  generator 


T 

K 

V 

m/s 

B 

T 

k 

_max 

q  g 

MW/m3 

600 

2000 

1 

0,5 

4,9 

600 

2000 

2 

0,5 

77,8 

600 

3000 

1 

0,5 

24,6 

600 

3000 

2 

0,5 

394,0 

600 

2000 

1 

0,2 

2,0 

600 

2000 

2 

0,2 

31,9 

300 

2000 

1 

0,5 

1,2 

300 

3000 

2 

0,5 

98,5 

Thus,  calculations  show  that  MHD  generator  with  non-equilibrium  conductivity  could  be 
realized  under  conditions  typical  for  HFV  flight,  with  self-sustained  mode  for  MHD  generator 
and  ionizer  being  subject  to  excessive  power  qg-  q,  constituting  considerable  part  of  MHD 
generator  power. 

In  closing  the  paragraph,  we  shall  consider  one  more  important  aspect  of  MHD  generator 
interaction  with  other  systems  on  board  HFV.  Since  such  generator  is  a  part  of  MPCE,  electric 
energy  production  is  not  a  single  function  it  should  perform.  In  line  with  AJAX  general  concept, 
MHD  generator  has  to  provide  for  efficient  control  of  MPCE  performance.  In  the  succeeding 
sections  an  important  MPCE  parameter  is  introduced  and  substantiated,  the  so  called  parameter 
of  MHD  interaction  Sv  =oB2L/pv,  where  L  is  MHD  generator  length.  As  it  will  be  shown  later 
on,  to  enhance  markedly  ramjet  performance,  parameter  of  MHD  interaction  has  to  meet 
inequality  Sv  >0.1.  As  per  this  inequality,  if  we  are  to  decrease  energy  supplied  to  flow  for  its 
ionization,  then  it  is  necessary  to  increase  the  value  of  magnetic  induction. 

Fig.6  depicts  region  of  parameters  B  and  q;  variation,  wherein  self-sustained  mode  for 
MHD  generator  with  non-equilibrium  conductivity  is  provided  at  the  same  time  with  required 
value  of  MHD  interaction  parameter.  Incident  curve  here  corresponds  to  the  permanent  value  of 
MHD  interaction  parameter  Sv=0.1.  Rising  curve  corresponds  to  critical  conditions  of  MHD 
generator  for  which  qg=qi.  Parameters  q;  and  B  corresponding  to  the  points  on  the  plane  in 
crosshatched  region  meet  simultaneously  self-sustained  mode  of  MHD  generator  with  non¬ 
equilibrium  conductivity  and  the  condition  of  MHD  interaction  Sv  >  0.1.  It  follows  from  the 
presented  figure  that  there  exists  optimal  value  of  energy  deposition  into  ionization  wherein 
existence  of  the  required  MPCE  operating  conditions  is  ensured  under  minimum  value  of 
magnetic  induction. 
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8.  Characteristics  of  MHD-controlled  inlet  in  a  wide  range  of  parameters  variation 
(altitude,  Mach  number,  inlet  geometry,  MHD  and  ionizer  parameters)  will  be 
calculated.  Requirements  for  magnetic  system  and  ionizer,  which  ensure  efficient 
control  of  inlet  at  off-design  conditions,  will  be  formulated 

As  shown  in  the  paper,  MPCE  investigation  has  been  conducted  for  the  selection  of  MHD 
channel  type  and  its  gasdynamic  profile  wherein  parameters  of  MHD  flow  ensure  maximum 
production  of  electric  energy  per  a  mass  unity  of  MHD  generator  under  preset  inlet  parameters 
and  restrictions.  When  calculating  mass,  only  components  of  the  channel  and  magnet  have  been 
taken  into  consideration.  The  investigation  was  conducted  by  the  example  of  two  types  of  flight 
vehicles  (FV)  -  a  full-scale  FV  and  a  small  model.  Various  alternatives  of  MHD  generator 
position  in  relation  to  combustion  chamber  have  been  considered.  In  case  of  the  full-scale  FV, 
the  generator  was  positioned  upstream  of  the  ramjet.  In  case  of  the  model  two  versions  of  MHD 
generator  arrangement  have  been  considered  -  upstream  of  the  ramjet  and  downstream  from  it. 

As  mentioned  earlier,  MHD  generator  positioned  upstream  of  the  ramjet  is  designed  for  FV  inlet 
control.  It  is  capable  of  electric  energy  production  which  can  be  used  on  board  FV  for  the 
solution  of  various  problems.  One  of  such  problems  is  electric  power  supply  to  ionizer  for  the 
generation  of  ionized  flow  aimed  at  provision  of  required  of  MHD  interaction  in  the  inlet.  Under 
high-velocity  flight  conditions  when  MHD  generator  produces  more  electric  power  than  is 
needed  for  ionizer  operation,  the  excess  energy  could  be  used  by  MHD  accelerator  -  a  special 
device  converting  electrical  power  into  kinetic  energy  of  ionized  gas  flow  at  the  nozzle  outlet. 
Excessive  energy  can  also  be  used  for  the  production  of  various  high-power  directed  beams  with 
the  aim  of  changing  properties  of  the  flow  around  the  flight  vehicle. 

The  investigation  conducted  to  optimize  MHD  generator  boiled  down  to  the  selection  of 
electrode  connection  diagram  for  MHD  channel  and  choosing  profile  of  gasdynamic  channel 
ensuring  maximum  of  the  produced  energy  per  a  mass  unit  at  the  specified  inlet  parameters  and 
restrictions.  As  a  first  approximation,  the  ratio  of  the  energy  being  produced  in  MHD  channel  to 
the  mass  of  magnetic  system  may  be  taken  as  specific  energy  of  MHD  generator  because  in  the 
configuration  of  MHD  generator  under  consideration  the  mass  of  magnetic  system  comprises 
around  80  %  of  the  total  mass  of  the  whole  energy  MHD  system.  For  the  purposes  of  this 
investigation  an  ultra-conductive  magnetic  system  has  been  selected  as  it  does  not  require  energy 
supply  under  the  operating  conditions.  Apart  from  that,  according  to  the  world  experience  on  the 
development  of  the  mobile  MHD  power  systems,  ultra-conductive  magnetic  systems  possess  the 
least  mass  as  compared  to  other  types  of  magnetic  systems.  It  was  also  presumed  that  making 
superconductive  magnetic  system  serviceable,  powering  and  cooling  inclusive,  is  exercised 
through  the  ground  auxiliary  complex  during  pre-flight  training. 

Two  types  of  MHD  generator  have  been  considered  -  with  the  Faraday  type  channel 
(with  continuous  electrodes)  and  the  Hall  MHD  channel  with  rectangular  section.  Selection  of 
these  channels  is  associated  with  the  optimum  effect  of  their  application  in  the  specified  region 
of  the  Hall  parameter  variation.  The  Faraday  MHD  channel  with  continuous  electrodes  differs 
advantageously  from  the  one  with  sectioned  electrodes  with  regard  to  load  conditions  as  in  the 
first  case  only  one  load  is  applied. 

The  following  restrictions,  common  to  all  cases,  have  been  imposed  on  geometrical 
dimensions:  the  channel’s  length  should  not  exceed  5m,  section  area  at  the  input  is  0.2  m2, 
coefficient  of  expansion  for  cross-section  area  should  not  exceed  2.  Electric  intensity  in  the 
MHD  channel  at  that  should  not  exceed  4  kV/m.  The  last  restriction  is  connected  with  the 
possibility  of  electrical  break-down  through  channel  walls. 

The  optimization  was  conducted  using  the  following  system  of  parameters: 

1)  Yi  —  distance  between  walls-electrodes  at  the  input  to  the  section  of  transformation  of 

MHD  channel; 

2)  Y2  -  distance  between  walls-electrodes  at  the  output  from  the  section  of  transformation  of 
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MHD  channel; 

3)  Zi  -  distance  between  walls-insulators  at  the  input  to  the  section  of  transformation  of 

MHD  channel; 

4)  Z2  -  distance  between  walls-insulators  at  the  output  from  the  section  of  transformation  of 

MHD  channel; 

5)  L  -  length  of  the  section  of  MHD  channel  transformation; 

6)  B  -  magnetic  field  induction  in  the  channel  of  the  section  of  MHD  transformation; 

7)  voltage  generated  by  MHD  generator  (for  the  Faraday  channel); 

8)  current  generated  by  MHD  generator  (diagonal  channel). 

The  procedure  of  optimization  run  as  follows: 

At  first,  using  established  restrictions  efficiency  function  was  being  formed.  Then,  based 
on  gasdynamic  calculations  and  using  initial  data,  energetic  parameters  and  dimensions  of  the 
working  section  of  MHD  channel  have  been  determined.  Calculations  were  carried  out 
numerically  to  a  quasi-one-dimensional  approximation.  It  was  based  on  the  approximation  of  the 
canon  flow  with  inhomogeneous  fluid  dynamics  parameters.  The  flow  showed  up  as  non-viscous 
non-heat-conductive  core  with  uniform  distribution  of  longitudinal  velocity  and  enthalpy 
surrounded  by  boundary  layers  on  the  walls  of  MHD  channel. 

Dimensions  of  working  section  of  MHD  channel  determine  the  volume  of  the  “warm” 
part  of  magnetic  system.  Starting  from  this  volume  and  quantity  of  magnetic  induction, 
dimensions  and  mass  of  magnetic  system  were  being  determined.  General  calculation  scheme 
was  built  up  in  such  a  way  that  the  parameters  of  a  magnetic  system  could  be  determined  also  for 
the  magnets  of  the  conventional  non-ultra-conductive  type. 

When  designing  a  full-scale  FY,  the  following  initial  data  and  limits  have  been  used: 

1.  The  flow  gas  is  the  air  of  the  following  composition  (in  volume  fractions):  N2  -  0.7826, 
02  -  0.2100,  C02  -  0.00030,  H2  -  0.00014,  inert  gases  -  0.00930. 

2.  Working  parameters  of  the  flow  at  MHD  generator  inlet:  static  pressure  -  1  and  3  atm; 
static  temperature  -  1500  K;  full  enthalpy  -  5,  10  and  15  MJ/kg;  electric  conductivity 
(usage  of  ionizer  is  stipulated)  -  1  cm/m. 

3.  Continuous  operation  time  -  2  h. 

The  following  parameters  were  being  optimized:  geometric  dimensions  (Yi,  Y2,  Zi,  Z2, 
L),  generator’s  power  Qg,  magnetic  induction  B.  Efficiency  function  was  being  defined  as 
maximum  of  the  ratio  generator’s  power  to  the  mass  of  magnet  (specific  power  of  generator  Qsp). 
Optimization  results  are  given  in  Tables  1  -  3  where  oce  is  electrode  plane  inclination  angle  in 
MHD  channel  of  diagonal  type;  h0  is  full  specific  enthalpy;  pi  and  P2  are  pressures  at  the  input 
and  output  of  the  section  of  MHD  channel  transformation,  respectively;  M  is  Mach  number,  m  is 
mass  of  magnet. 

It  follows  from  the  data  presented  that  optimal  dimensions  of  gasdynamic  channel  are 
closely  allied  practically  for  all  combinations  of  initial  data  for  MHD  channel  of  this  type. 

MHD  channel  dimensions  have  been  obtained  proceeding  from  the  fact  that  under 
conditions  of  low  electric  conductivity  of  the  working  gas  (1  cm/m)  and,  consequently,  low 
value  of  MHD  interaction  quantity,  MHD  channel  tends  to  acquire  such  geometrical  dimensions 
within  specified  restrictions  that  would  ensure  its  maximum  volume.  It  has  to  do  with  MHD 
generator  power  being  proportional  to  the  volume  of  working  section  of  MHD  channel. 

From  the  results  obtained  follows  that  apart  from  the  point  (ho  =  10  MJ/kg,  pi  =  1.0  atm.) 
for  the  Faraday  channel  where  other,  less  optimal  local  minimums  have  also  been  found,  optimal 
value  of  Ng  does  not  depend  on  the  inlet  pressure  of  channel  of  MHD  transformation  section. 
This  value  depends  only  on  full  enthalpy  of  the  working  flow.  For  all  parameter  combinations, 
save  the  above,  magnetic  induction  at  optimum  conditions  makes  0.55  T.  It  has  to  do  with 
magnetic  system  mass  being  directly  proportional  to  the  induction  value  of  magnetic  field. 

The  same  results  have  been  obtained  for  MHD  channel  of  the  Hall  type  but  here  optimal 
induction  value  of  magnetic  field  depends  on  the  inlet  pressure  of  channel  of  MHD 
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transformation  section.  Really,  with  pressure  increase  electrons  mobility  drops  and  to  obtain 
optimal  value  of  all  parameter  it  is  necessary  to  increase  induction  value  of  magnetic  field. 


Table  1.  Optimal  values  of  geometrical  parameters  for  a  full-scale  FV 


MHD  channel  type 

L,  m 

Yi,m 

Y2,m 

Zh  m 

Z2,  m 

a,E,  degree 

Faraday 

5.000 

0.479 

0.710 

0.400 

0.560 

Hall 

5.000 

0.591 

0.790 

0.340 

0.500 

10 

Table  2.  General  results  of  parameter  optimization  for  a  full-scale  FY  with  Faraday 
MHD  channel 


h0 

Pi 

P2 

M 

Qg 

m 

QsP 

B 

MJ/kg 

atm 

atm 

MW 

kg 

kW/kg 

T 

5.0 

1.0 

0.4 

3.4 

0.87 

300 

2.90 

0.55 

3.0 

1.2 

3.36 

0.87 

300 

2.90 

0.55 

10.0 

1.0 

0.4 

5.20 

5.12 

1800 

2.84 

1.35 

3.0 

1.2 

5.16 

1.99 

300 

6.63 

0.55 

15.0 

1.0 

0.4 

6.50 

3.27 

300 

10.9 

0.57 

3.0 

1.2 

6.50 

3.10 

300 

10.3 

0.55 

Table  3.  General  results  of  parameter  optimization  for  a  full  -scale  FY  with  Hall 
Hall  MHD  channel 


ho 

Pi 

P2 

M 

Qg 

m 

Qsp 

B 

MJ/kg 

atm 

atm 

MW 

kg 

kW/kg 

T 

5.0 

1.0 

0.4 

3.4 

1.50 

600 

2.50 

0.764 

3.0 

1.2 

3.4 

4.19 

1600 

2.60 

1.27 

10.0 

1.0 

0.4 

5.23 

3.06 

500 

6.12 

3.0 

1.2 

5.16 

10.30 

1700 

6.05 

15.0 

1.0 

0.4 

6.58 

5.51 

600 

9.18 

0.764 

3.0 

1.2 

6.49 

16.20 

1700 

9.50 

1.300 

Optimal  conditions  of  pressure  variation  mode  in  MHD  channel  have  been  found  for  all 
cases  under  investigation. 

And  now  let  us  switch  to  the  description  of  MHD  generator  optimization  results  for  a 
small  FV  model.  Let  us  consider  a  case  where  MHD  generator  is  arranged  along  the  flow 
upstream  of  combustion  chamber. 
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Like  in  the  previous  case,  the  flow  gas  is  the  air  of  following  composition  (in  volume 
fractions):  N2  -  0.7826,  O2  -  0.2100,  CO2  -  0.00030,  H2  -  0.00014,  inert  gases  -  0.00930 

Working  parameters  of  the  flow  at  MHD  generator  inlet  have  the  following  values:  static 
pressure  is  0.5  atm,  static  temperature  -1500  K,  flow  velocity  at  the  inlet  to  the  section  of  MHD 
transformation  channel  -  3000  m/s;  specific  electrical  conductivity  (usage  of  ionizer  is 
stipulated)  -  lcm/m.  An  additional  restriction  has  been  established  -  magnetic  field  induction 
should  not  exceed  2  T. 

The  same  as  with  a  full-scale  FV,  subjected  to  optimization  were  geometrical  dimensions 
(Yi,  Y2,  Zi,  Z2,  L),  power  of  the  generator  Qg,  magnetic  induction  B.  Specific  power  Qsp  was  the 
efficiency  function. 

The  results  of  optimization  are  presented  in  Tables  4-5.  Accepted  designations  now  and 
hereinafter  are  Ai  and  A2  -  section  areas  at  the  inlet  and  outlet  to  and  from  the  channel  of  MHD 
generator  transformation,  respectively. 

The  same  as  in  the  previous  case,  it  follows  from  the  results  presented  that  under  the 
established  geometrical  limits  MHD  channel  seeks  to  occupy  the  fullest  volume  due  to  the  low 
value  of  flow  electric  conductivity  (1  cm  m"1). 

Table  4.  Optimal  values  of  MHD  channel  geometrical  parameters  for  a  small  FV 


model  at  MHD  generator  location  upstream  of  combustor 


MHD  channel 
type 

L,  m 

Yum 

Y2,  m 

Zi,m 

Z2,  m, 

Ai/A2 

Faraday 

5.000 

0.50 

0.70 

0.40 

0.57 

2.00 

Hall 

5.000 

0.50 

0.45 

0.40 

0.84 

1.87 

Table  5.  General  results  of  optimization  of  Faraday  and  Hall  MHD  channels  for 
a  small  FV  model  at  MHD  generator  location  upstream  of  combustor 

MHD  channel  pi 

type  atm 

p2  M 

atm 

Qg 

MW 

m 

kg 

Qsp  B 

kW/kg  T 

Faraday  0.50 

0.23  3.7 

6.57 

5594 

1.10  2.0 

Hall  0.50 

0.25  3.7 

8.56 

6700 

1.28  2.0 

Optimal  set  of  parameters  for  both  MHD  channel  types  has  been  obtained  for  the  lowest 
value  of  magnetic  induction  specified  by  the  above  restriction.  As  in  the  case  with  a  full-scale 
FV,  it  can  be  explained  by  the  quest  for  the  least  magnet  mass.  But  as  in  this  case  magnetic 
induction  was  limited  by  the  low  value,  optimal  values  of  generator’s  power  and  magnet  mass 
turned  to  be  higher  and  specific  power  of  generator  lower  than  for  a  full-scale  model.  Besides, 
optimal  values  of  generator’  power  and  magnet’  mass  were  found  to  be  higher  for  the  Hall 
channel  than  those  for  the  Faraday  channel. 

Optimal  values  of  pressure  parameters  in  MHD  channel  have  also  been  determined. 

And  now  let  us  turn  to  the  description  of  MHD  generator  optimization  results  for  a  small 
FV  model  in  the  case  of  MHD  generator  arrangement  along  the  flow  downstream  of  combustor. 

Three  modes  of  operation  have  been  considered  lasting  20,  60,  and  120  sec  respectively. 

The  flow  gas  comprised  kerosene  combustion  yields  (the  ratio  kerosene-air  being  0.1). 
The  working  flow  parameters  at  MHD  generator  inlet  had  the  following  values:  static  pressure  - 
0.5  atm,  static  temperature  -  2500  K;  inlet  flow  velocity  at  the  section  of  MHD  transformation 
channel  -  3000  m/s. 

The  following  additional  restrictions  have  been  imposed:  magnetic  induction  should  be 
no  less  than  2  T  and  generator’  power  -  no  less  than  10  MW.  Subjected  to  optimization  were 
geometrical  dimensions  (F/,  Y2,  Zu  Z2,  L),  power  of  the  generator  Qg  and  magnetic  induction  B. 
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Thermal  ionization  of  the  medium  has  been  taken  into  consideration  and  to  enhance  ionization 
degree  K-Na  eutectic  was  introduced  into  the  flow.  Consumption  rate  of  the  eutectic  was  also  the 
optimization  parameter.  Minimum  sum  of  magnet  mass  m  and  mass  of  easily-ionized  additive  ms 
was  viewed  in  this  case  as  efficiency  function. 

Optimization  results  are  presented  in  Tables  6-8.  Here,  eis  percentage  of  additive  in  the 

flow. 

From  the  results  presented  it  follows  that  MHD  channels  in  this  particular  case  are  shorter 
and  more  converging  which  may  be  attributed  to  higher  conductivity  of  the  working  flow  (~  25 
cm/m)  and  higher  value  of  MHD  interaction  parameter,  respectively. 

The  optimal  set  of  parameters  for  both  Faraday  and  Hall  types  of  channels  was  obtained 
at  minimum  permissible  power  of  the  generator  (10  MW)  and  induction  value  of  magnetic  field 
(2  T).  The  value  of  specific  power  of  the  generator  is  decreased  with  operation  time  increase  that 
can  be  explained  by  quantitative  increase  in  the  mass  of  the  additive. 

Table  6.  Optimal  values  of  MHD  channel  geometrical  parameters  for  a  small  FV 


model  in  case  of  MHD  generator  location  downstream  of  combustor 


MHD  channel 

L, 

Yi, 

y2, 

Zi, 

z2, 

a,/a2 

type 

m 

m 

m 

m 

m 

Faraday 

0.90 

0.500 

0.308 

0.400 

0.202 

0.31 

Hall 

0.46 

0.500 

0.570 

0.400 

0.220 

0.65 

Table  7.  General  optimization  results  of  the  Faraday  MHD  channel  for  a  small  FV 
model  in  case  of  MHD  generator  location  downstream  of  combustor 


top 

Pi 

P2 

M 

Qg 

m 

m  +  ms 

Qsp 

B 

8 

s 

atm 

atm 

MW 

kg 

kg 

kW/kg 

T 

% 

20 

0.5 

0.3 

3.35 

10 

1894 

1927 

5.20 

1.97 

4.18 

60 

0.5 

0.3 

3.35 

10 

1894 

1994 

5.024 

1.97 

4.18 

120 

0.3 

0.3 

3.35 

10 

1894 

2095 

4.80 

1.97 

4.18 

Table  8.  General  optimization  results  of  the  Hall  MHD  channel  for  a  small  FV 


model  in  case  of  MHD  generator  location  downstream  of  combustor 


top 

Pi 

P2 

M 

Qg 

m 

m  +  m 

Qsp 

B 

8 

S 

atm 

atm 

MW 

kg 

kg 

kW/kg 

T 

% 

20 

0.5 

0.98 

3.35 

10 

1758 

1784 

5.6 

2.07 

3.2 

60 

0.5 

0.97 

3.35 

10 

1762 

1840 

5.4 

2.07 

3.2 

120 

0.5 

0.97 

3.35 

10 

1753 

1908 

5.2 

2.07 

3.2 

Some  higher  value  of  specific  power  was  obtained  for  the  Hall  MHD  generator. 


Thus,  optimized  (as  per  Ajax  conception)  MHD  generators  designed  for  flight  vehicles 
have  specific  power  around  2.9  -  1 1  kW/kg  at  full  power  0.87  -  5.12  MW  for  MHD  channel  of 
the  Faraday  type.  For  the  Hall  type  generators  specific  power  is  within  2.5  -  9.5  kW/kg  range 
and  full  power  is  within  1.5  -  16.2  MW  range.  But  optimal  conditions  of  Hall  generator  are 
realized  at  greater  mass  of  the  magnet. 


70 


Conclusion 


The  investigation  on  feasibility  of  MHD  system  principles  application  in  hypersonic 
technologies  for  control  of  a  HFV  aerodynamic  performance  has  proved  reality  of  their 
realization  when  creating,  in  particular,  an  engine  of  the  latest  design  -  MPCE,  as  well  as 
developing  a  number  of  engineering  solutions,  inclusive  of  the  following: 

•  conception  of  MHD  systems  application  for  the  air  flow  deceleration  and  scramjet 
efficiency  enhancement; 

•  conception  of  MHD  generator  application  for  control  of  mass  flow  rate  and  the  flow 
structure  in  the  scramjet  air  inlet; 

•  conception  and  mathematical  model  of  MHD  generator  on  cold  air  flow  with  non¬ 
equilibrium  ionization. 

The  elaboration  of  the  conception  of  MHD  systems  application  for  deceleration  and 
compression  of  the  air  flow  as  well  as  for  the  scramjet  efficiency  enhancement  has  been  carried 
out  with  the  analysis  of  MPCE  system  consisting  of  the  air  inlet,  an  internal  MHD  generator  of 
the  Faraday  type,  combustion  chamber,  MHD  accelerator  and  a  nozzle.  General  expressions  have 
been  evaluated  for  specific  impulse  which  is  a  function  of  many  MPCE  parameters.  The  area  of 
their  variation  has  been  thoroughly  studied  wherein  MHD  system  ensures  increase  of  the 
propulsion  specific  impulse.  Its  dependence  on  load  factor  of  MHD  generator  has  a  maximum 
corresponding  to  the  inlet  with  a  minimum  value  of  the  angle  of  flow  turn,  ft  was  shown  that  a 
prerequisite  to  MPCE  impulse  enhancement  is  pressure  increase. 
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